THE 


PROCEEDINGS 
OF 


fie EtYOICAL. SOCIETY 


Section B 
Voi. 65, PART6 1 June 1952 No. 390B 
CONTENTS 
PAGE 
Dr. P. T. LanpsBerc. Further Results in the General peay of Barrier Layer 
Rectifiers : : : : f ‘ é 397 
Dr. R. Cooper. Breakdown in Selenium Rectifiers . . 409 
Dr. G. F. J. Gartick and Mr. G. T. Wricur. Chinlectorcies of Scintillation 
Counters f J é 3 ; f : 415 
Dr. L. Jacos. The Current in the Electron Immersion Obiesave 421 
Dr. D. G. Avery. An kaa Method for Measurements of tte? cee 
by Reflection 425 
Mr. C.G. WYNNE. Primary Ana and Cokes Change! . 429 
Dr. H. D. Parsrook and Dr. E. G. RICHARDSON. ple aes of Ultrasonic Waves 
in Vapours near the Critical Point : : : : . 437 
Prof. E. N. pA C. ANDRADE and Dr. R. F. Y. Bye, The Influence of Electro- 
lytes on the Mechanical Properties of Certain Metal Single Crystals . 445 
Letters to the Editor : 
Mr. G. D. Apam and Dr. K. J. StanpLey. Ferromagnetic Resonance in 
Manganese Antimonide 454 
Mr. E. W. Lez. Low Remanence and the Temperature Variation of Perme- 
ability of Silicon Iron Alloys. 455 
Mr. J. M. Natsu and Mr. E. R. Wess. Strain Patterns in Toughened Ciscoe . 457 
Mr. J. Mazur. On the Sampling of Water Droplets in Natural Clouds and in 
Radiation Fogs . i } : y : » 457 
Mr. P. F. Litrve and Dr. A. VON ENGEL. The Dispersion of Electron Beams 
in Gases q 459 
Mr. F. F. ROBERTS. Microwave F araday Rotation in Liquid Oxygen 460 
Dr. R. Street, Dr. J. C. Woo.iey and Mr. P. B. Smirn. The Influence of 
Heat Treatment on Magnetic Viscosity in Permanent Magnet Alloys 461 
Contents for Section A 462 
Abstracts for Section A 463 


Price to non-members 10s. net, by post 9d. extra. Annual subscription: £5 5s. 


Composite subscription for both Sections A and B: £9 9s. 


Published by 
THE PHYSICAL SOCIETY 
1 Lowther Gardens, Prince Consort Road, London S.W.7 


THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


PROCEEDINGS OF THE PHYSICAL SOCIETY 


The Proceedings is now published monthly in two Sections. 


ADVISORY BOARD 


Chairman: The President of the Physical Society (R. WHIDDINGTON, C.B.E., M.A., D.Sc., F.R.S.) 


E. N. pa C. ANpRADE, Ph.D., D.Sc., F.R.S. 

Sir Epwarp APPLETON, G.B.E., K.C.B., 
D.Sc., F.R.S. 

P. M. S. BLAcKeTT, M.A., F.R.S. 

Sir-LAWRENCE Bracc, O.B.E., M.C., M.A., 
Sc.D., D.Sc., F.R.S. 

Sir James Cuapwick, D.Sc., Ph.D., F.R.S. 

S. CuapmMan, M.A., D.Sc., F.R.S. 

ene oF OxrorpD, M.A., Ph.D., 

Sir JoHN CockcrorT, C.B.E., M.A., Ph.D., 
F.R.S. 


Sir Cuartes Darwin, K.B.E., M.C., M.A., 
Sc.D., F.R.S. 

N. FeaTuer, Ph.D., F.R.S. 

G. I. Fincu, M.B.E., D.Sc., F.R.S. 

D. R. Hartree, M.A., Ph.D., F.R.S. 

N. F. Mort, M.A., D.Sc., F.R.S. 

M. L. OLIPHANT, Ph.D., D.Sc., F.R.S. 

F. E. Srwon, C.B.E., M.A., D.Phil., F.R.S. 

T. SmitH, M.A., F.R.S. 

Sir GEORGE THOMSON, M.A., D. Sc., F.R.S. 


Papers for publication in the Proceedings should be addressed to the Hon. Papers Secretary, 
Dr. H. H. Hopkins, at the Office of the Physical Society, 1 Lowther Gardens, Prince 
Consort Road, London S.W.7. Telephone : KENsington 0048, 0049. 


Detailed Instructions to Authors can be obtained from the Secretary—Editor. 


BULLETIN ANALYTIQUE 
Publication of the Centre National de la Recherche Scientifique, France 

The Bulletin Analytique is an abstracting journal which appears in three parts, Part 1 
covering scientific and technical papers in the mathematical, chemical and physical sciences and their 
applications, Part 2 the biological sciences and Part 3 philosophy. 

The Bulletin, which started on a modest scale in 1940 with an average of 10,000 abstracts per 
part, now averages 35 to.45,000 abstracts per part. The abstracts summarize ‘briefly papers in 
scientific and technical periodicals received in Paris from all over the world and cover the majority 
of the more important journals in the world scientific press. The scope of the Bulletin is constantly | 
: being enlarged to include a wider selection of periodicals. 

The Bulletin thus provides a valuable reference book both for the laboratory and for the individual 
research worker who wishes to keep in touch with advances in subjects bordering on his own. 

A specially interesting feature of the Bulletin is the microfilm service. A microfilm is made of 
each article as it is abstracted and negative microfilm copies or prints from microfilm can be purchased 
from the editors. 

The subscription rates per annum for Great Britain are 4,000 frs. (£4) each for Parts 1 and 2, and 
74 ae frs. (£2) for Part 3. Subscriptions can also be taken out to individual sections of the Bulletin 
as follows : 


Pure and Applied Mathematics—Mathematics—Mechanics 

Astronomy—aAstrophysics—Geophysics . 

General Physics—‘Thermodynamies-—Heat—Optics—Elec- 
tricity and Magnetism 

Atomic Physics—Structure of Matter 

General Chemistry—Physical Chemistry 

Inorganic Chemistry — Organic Chemitey => “Applied 
Chemistry—Metallurgy : 

Engineering Sciences . . 

Mineralogy-Petrography—Geology—Palacontology 

Biochemistry—Biophysics—Pharmacology 

Microbiology—Virus and Phages : 

Animal Biology—Genetics—Plant Biology ae 

Agriculture—Nutrition and the Food Industries . . 


Subscriptions can be paid directly to the editors : Centre Nees ae la Recherche Scientifique, 
18, rue Pierre-Curie, Paris 5¢me (Compte-chéque-postal 2,500-42, Paris), or through Messrs. H. K. 
Lewis & Co. Ltd. , 136 Gower Street, London W.C.1. 


THE PROCEEDINGS OF THE PHYSICAL SOCIETY i 


THE PHYSICAL SOCIETY 


MEMBERSHIP 

Membership of the Society is open to all who are interested in Physics: 

FELLOwsHIP. A candidate for election to Fellowship must as a rule be recommended 
by three Fellows, to two of whom he is known personally. Fellows may attend all meetings 
of the Society, are entitled to receive Publications 1 (either Section A or Section B), 4 and 5 
below, and may obtain the other publications at much reduced rates. 

STUDENT MEMBERSHIP. A candidate for election to Student Membership must be 
between 18 and 26 years of age and must be recommended from personal knowledge by 
a Fellow. Student Members may attend all meetings of the Society, are entitled to receive 
Publications 1 (either Section A or Section B) and 4, and may obtain the other publications 
at much reduced rates. 

Books and periodicals may be read in the Society’s Library, and a limited number of 
books may be borrowed by Fellows and Student Members on application to the Honorary 
Librarian. 

Fellows and Student Members may become members of the Colour Group, the Optical 
Group, the Low Temperature Group and the Acoustics Group (specialist Groups formed in 
the Society) without payment of additional annual subscription. 


PUBLICATIONS 


1. The Proceedings of the Physical Society, published monthly in two Sections, contains 
original papers, lectures by specialists, reports of discussions and of demonstrations, and 
book reviews. Section A contains papers mainly on atomic and sub-atomic subjects; Section B 
contains papers on macroscopic physics. 

2. ‘Reports on Progress in Physics, published annually, is a comprehensive review by 
qualified physicists. 

3. The Handbook of the Physical Society’s Annual Exhibition of Scientific Instruments 
and Apparatus. ‘This Exhibition is recognized as the most important function of its kind, 
and the Handbook is a valuable book of reference. 

4. The Bulletin, issued at frequent intervals during the session, informs members of 
programmes of future meetings and of the business of the Society generally. 

5. Physics Abstracts (Science Abstracts A), published monthly in association with the 
Institution of Electrical Engineers, covers the whole field of contemporary physical research. 

6. Electrical Engineering Abstracts (Science Abstracts B), published monthly in association 
with the Institution of Electrical Engineers, covers the whole field of contemporary research 
in electrical engineering. 

7. Special Publications, critical monographs and reports on special subjects prepared 
by experts or committees, are issued from time to time. 


MEETINGS 


At approximately monthly intervals throughout each annual session, meetings are 
held for the reading and discussion of papers, for lectures, and for experimental demonstra- 
tions. Special lectures include: the Guthrie Lecture, in memory of the founder of the 
Society, given annually by a physicist of international reputation; the Thomas Young 
Oration, given biennially on an optical subject; the Charles Chree Address, given biennially 
on Geomagnetism, Atmospheric Electricity, or a cognate subject; and the biennial Ruther- 
ford Memorial Lecture. Meetings are generally held each year at provincial centres, and 
from time to time meetings are arranged jointly with other Societies for the discussion of 


subjects of common interest. 
Each of the four specialist Groups holds about five meetings in each session. 


SUBSCRIPTIONS 


Fellows pay an Entrance Fee of £1 1s. and an Annual Subscription of £3 3s. Student 
Members pay only an Annual Subscription of 15s. Second Section of Proceedings 30s. 
No entrance fee is payable by a Student Member on transfer to Fellowship. 


Further information may be obtained from the Secretary-Edttor 
at the Offices of the Society : 


1 LowrnHerR GARDENS, PRINCE Consort Roap, Lonpon S.W. 7 
Telephone: KENsington 0048, 0049 


a 


il THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


atomic 


Fair Copying Declaration 


The Physical Society is one of the subscribe 
to the Fair Copying Declaration. Accordin 
to this it is permissible for single reproduction 
of articles published in its journal to be mad 
for the purpose of private study or research 
without specific application to the ae | 

| 


scientists 
AA 


The Official Journal of the ATOMIC SCIENTISTS’ ASSOC. 


@ The Atomic Scientists’ Association is an association of 
scientists whose work has given them special know- 
ledge of the consequences for the world of the use and 
misuse of atomic energy. To make known the true 
facts about atomic energy and its implications it pub- 
lishes every two months the Atomic Scientists’ News. 


The following regulations cover the agreement 
1. The recipient of the copy is given notice 
that he is liable for infringement of copy’ 
right by misuse of the copy, and that it i 
illegal to use the copy for any furthe 


CONTENTS for No. 5. Vol. I. reproduction. 
Symposium on Civil Defence : 


2. The organization making and furnishing tha 


Assessment of Atomic Casualties. By E. C. Allen. 

Radioactivity and Civil Defence. By Dr. W. G. Marley. . : 

The Biological Component. By D. G. Arnott. copy does so without profit to itself. 
The Casualty Service. By Sir Claude Frankau. 5 


Administration and Organisation. By S. L. Harford. 


Prospects for the United Nations Disarmament 3. Proper acknowledgment Se. to the 


Commission. By Robert W. Frase. publication from which the copy is made. 
Review of the Month. By F. R. N. Nabarro 
Book Reviews. Prof. Sir George Thomson, F.R.S. 
Dr. G. Stephenson. 4. Not more than one copy of any one excerpt: 
PRICE 6/- per part Ann. Sub. £1-12-6 post free shall be furnished to any one person. 


Printed and published by 
TAYLOR & FRANCIS, LTD., Red Lion Ct., Fleet St., E.C.4 


PHYSICAL SOCIETY. SPECIALIST GROUK 


OPTICAL GROUP 

The Physical Society Optical Group exists to foster interest in and development of all branches « 
optical science. ‘To this end, among other activities, it holds meetings about five times a year t 
discuss subjects covering all aspects of the theory and practice of optics, according to the papers offerec 
COLOUR GROUP 

The Physical Society Colour Group exists to provide an opportunity for the very varied types « 
worker engaged on colour problems to meet and to discuss the scientific and technical aspects of the: 
work. Five or six meetings for lectures and discussions are normally held each year, and reprints ¢ 
papers are circulated to members when available. A certain amount of committee work is undertakes 
and reports on Defective Colour Vision (1946) and on Colour Terminology (1948) have already bee 
published. 
LOW TEMPERATURE GROUP 

The Low Temperature Group was formed to provide an opportunity for the various groups « 
people concerned with low temperatures—physicists, chemists, engineers, etc.—to meet and becor 
familiar with each other’s problems. The group seeks to encourage investigations in the low temperatu1 
field and to assist in the correlation and publication of data. 
ACOUSTICS GROUP 

The Acoustics Group was formed to meet the long-felt need for a focus of acoustical studies i 
Great Britain. ‘The scope includes the physiological, architectural, psychological, and musical aspect 
of acoustics as well as the fundamental physical studies on intensity, transmission and absorption « 
sound. ‘Ihe Group achieves its object by holding discussion meetings, by the circulation of reprint 
and by arranging symposia on selected acoustical topics. 


Further information may be obtained from the Offices of the Society : 
1 LOwTHER GARDENS, PRINCE ConsorT RoaD, LONDON S.W.7. 


THE PROCEEDINGS OF THE PHYSICAL SOCIETY iti 


PROCEEDINGS OF THE 
PHYSICAL SOCIETY 


ADVERTISEMENT RATES 


The Proceedings are divided into two 
parts, Aand B. The charge for insertion 
is £18 for a full page in either Section A 
or Section B, £30 for a full page for 
insertion of the same advertisement in 
both Sections. The corresponding 
charges for part pages are: 


$page £9 5 0 7,15 2107 0 
tpage £415 0 £8 505.0 
$ page £210 0 £4 5 0 


Discount is 20% for a series of six 
similar insertions and 10% for a series 
of three. 


The printed area of the page is 
84” x54”, and the screen number is 100. 


Copy should be received at the Offices 
of the Physical Society six weeks before 
the date of publication of the Proceedings. 


VACUUM 


A NEW SCIENTIFIC QUARTERLY 


©PROFOUNDLY important developments 
in the field of industrial and scientific effort follow 
hard on advances in high vacuum techniques. 
“Vacuum”, now in its second year, is therefore 
published as a unique quarterly to report these 
developments, the basic techniques and equipment 
for which are of interest in widely separated fields 
having no other common literature. 


The authoritative articles and news columns, 
together with the international literature abstracts 
extensively classified and specially printed for detached 
filing, provide indispensable information for scientists 
and industrialists who wish to keep abreast of the 
times in an essential rapidly expanding technology. 


Subscribers in 33 Countries 
OO 
All communications to: —THE EDITOR, VACUUM, 


WORSLEY BRIDGE ROAD, LONDON, S.E.26 
Telephone: Sydenham 4444 


W. EDWARDS & CO. (LONDON) LTD., ENGLAND 


WHESPHYSICAT, SOCIETY, 
VOLUME XV of the REPORTS ON PROGRESS IN PHYSICS 


A comprehensive annual review by specialist authors. The contents are as follows: 


L. Kettner. The Hydrogen Bond. 
K. I. Mayne. Mass Spectrometry. 


A. E. Taytor. 


Range—Energy Relations. 


R. S. Curistran. The Analysis of High Energy Neutron—Proton and 
Proton—Proton Scattering Data. 


G. K. BatcueLtor. ‘Turbulent Motion. 


A. FarrweaTuHer, F. F. Roperts and A. J. E. WELCH. 


Ferrites. 


. V. Raynor. The Band Structure of Metals. 


G 
D. K. C. MacDonatp and K. SARGINSON. 
R 


Galvanomagnetic Effects in 
Conductors. 


. Komprner. Travelling-Wave Tubes. 


Price 50s.; by post 51s. 6s. 
Members 27s. 6d.; by post 29s. 


Further information can be obtained from 
THE PHYSICAL SOCIETY 


1 Lowther Gardens, Prince Consort Road, London S.W.7 
Neer 


iv THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


SegaUAN AUNT AR emu NN ea ENON N MUNNAR ANA eH 


THE SCIENTIFIC WOk 


of 
RENE DESCARTES | 


1596 — 1650 


OESIENW 
ENGINEERS/ 


Seamless one-piece 


METAL BELLOWS 


By J. F. SCOTT, B.a., M.sc., PELD 


This book puts the chief mathematical and ph 
discoveries of Descartes in an accessible form, an 
an outstanding gap upon the shelf devoted to 


2.» container Vcaled: HYDRAULICALLY FORMED history of philosophy and science. 
hermetica BY A PROCESS UNIOUE There is to be found in this volume the consider 
gland IN THIS COUNTRY contribution that Descartes made to the phy, 


sciences, which involved much accurate wor 
geometrical optics and its bearing upon the prac 
problem of fashioning lenses, as also the de' 
problems of light and sight and colour. The cai 
treatment that Dr. Scott has accorded to. the wor! 
Descartes is welcome. The book is well worth reac 
and will be an asset to all libraries. This publicatic 
recommended and approved by the: Publication F 
Committee of the University of London. 


FOR Automatic coolant regulation. Movement for pressure 
change. Packless gland to seal spindle in high vacua. Reservoir 
to accept liquid expansion. Dashpot or delay device. Barometric 
measurement or control. Pressurised couplings where vibration 
or movement is present. Dust seal to prevent ingress of dirt. 
Pressure reducing valves. Hydraulic transmission. Distance 
thermostatic control. Low torque flexible coupling. Pressure 
sealed rocking movement. Pressurised rotating shaft seals. 
Aircraft pressurised cabin control. Refrigeration expansion 
valves. Thermostatic Steam Traps. Pressure amplifiers. Dif- 
ferential pressure measurements. Thermostatic operation of 
louvre or damper. Write for List No, V. 800-1. 


Drayton METAL BELLOWS 


Drayton Regulator and Instrument Co. Ltd., West Drayton, Middlesex 


212 pages, 7” x 10”, amply illustrated. 
PRICE £1-0-0 net To be published July 


Printed and Published by 
TAYLOR & FRANCIS LTD., Red Lion Court, Fleet St., London, 


Ee Aaulanlanian\anlselanianian(enanjanianianialan se ieuisuianianjanisalanianjeniaoiey 1 


HANDBOOK REPORTS ON 
OF THE PROGRESS IN PHYSICS 


PHYSICAL SOCIETY’S 
36th EXHIBITION Volume XIV (1951) 


OF 

* pacers Aas pees bea Absorption S Dectroscopa 
T by W. C. Price. Anomalous Fine Structure of Hydregen 
SCIENTIFIC IN STRUMENTS and Singly Ionized Helium, by W.E.Lamb, Jr. New 
Techniques in Optical Interferometry, by H. Kuhn. 
AND APPARATUS The Diffraction Theory of Aberrations, by E. Wolf. The 
Spectrum of the Airglow and the Aurora, by A. B. 
Meinel. The Microphysics of Clouds, by B. J. Mason 
and F. H. Ludlam. Angular Correlations in Nuclear 
1 952 Reactions, by M. Deutsch. The Nuclear Interactions 
of Cosmic Rays, by G. D. Rochester and W. V. G. 
Rosser. Nuclear Effects in Atomic Spectra, by E. W. 
Foster. A Critical Survey of Ionospheric Temperatures, 
by N. C. Gerson. Some Applications of Nuclear 

Physics in Medicine, by W. V. Mayneord. 


; by post 7s. 3d. 
Price £2 10s., postage and packing 1s. 6d. 


Orders, with remittances, to Orders, with remittances, should be sent to the publishers : 
THE PHYSICAL SOCIETY 


THE PHYSICAL SOCIETY 
1 Lowther Gardens, Prince Consort Road 


1 Lowther Gardens, Prince Consort Road, 
London S.W.7 London S.W.7 


THE PROCEEDINGS OF 
THE PHYSICAL SOCIETY 


Section B 


VoL. 65, Part 6 1 June 1952 No. 390 B 


Further Results in the General Theory of 
Barrier Layer Rectifiers 


By Po TS DAN DSBERG 
Department of Natural Philosophy, The University, Aberdeen 


MS. received 27th September 1951, and in final form 7th January 1952 


ABSTRACT. It is pointed out that the temperature variation Vp=—a—bT of the diffusion 
potential, as deduced from the forward characteristic of rectifiers, is in broad agreement 
with theory. Hence formulae are derived for the temperature variation of the zero resistance 
of a rectifier having an arbitrary distribution of impurity centres, and subject to either diode 
or diffusion theory. Interpretation of experiments with the aid of these formulae leads to 
the temperature dependence of (a) the effective mass of current carriers in germanium 
rectifiers, and (6) of the mopility of current carriers in selenium rectifiers. The latter 
curve has the same form as that obtained from Hall effect measurement, and it is pointed out 
that it is of broadly the same shape as the corresponding curves for high resistivity samples of 
silicon and germanium. It is shown that two d.c. characteristics of the same rectifier, but 
corresponding to different temperatures, can intersect. Formulae for the applied voltage and 
the temperature rise which lead to the onset of thermal instability are deduced for all 
impurity centre distributions, and for diode and diffusion theory. A general theory of current 
creep in rectifiers is outlined, and good agreement obtained with experiments on selenium 


rectifiers. 


§1. INTRODUCTION 


HE general theory of barrier layer rectifiers has been introduced so that 

a single analysis can cover all likely impurity centre distributions in the 

semiconductor. ‘The theory can be applied to experimental material of 
considerable variety. ‘There seems to be little point, however, in attempting 
an exhaustive enumeration of possible applications, since these should become 
evident as new experimental material becomes available. It may perhaps be of 
interest, however, to discuss matters (described by the headings of the six 
succeeding sections) which have recently been the subject of enquiry from the 
point of view of this general theory, and in this way to carry them somewhat 
further. 
The notation introduced previously (Landsberg 1951 b, c; these papers will 

be referred to as I and II respectively) will be adopted. 
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§2. THE TEMPERATURE DEPENDENCE OF BARRIER HEIGHT 
AND RECTIFIER RESISTANCE 


Consider a barrier layer having an arbitrary distribution m(x) of fixed space : 


charge. The potential energy of a current carrier in the barrier has its maximum 
at a plane whose distance 6 from the contact is given by (II, eqn. (6)) 
1 ra 4e _., 
p= 16 | n(x) dx— GVA). ioe SG 
The coordinate x =d marks the end of the barrier and « is the dielectric constant 
of the bulk semiconductor. The current density 7 flowing through the barrier 


when the voltage U/e is applied is given in diffusion theory by (II, eqn. (11)) 
J =Jovr?er(1 em euinty. where y =e?/2ebkT 


and ae Ghe aves (#2) exp 2 ie eee (2) 


for reverse and very small forward potentials. For the forward direction U is 
positive and 7 is negative since the current flows in the direction of the negative — 


x.axis. For the reverse direction U is negative and 7 is positive. 


If by be the value of 6 when U=0, then the corresponding maximum | 


potential energy in the barrier is (II, eqn. (7)) 


V(b) =Vpaej2eby. oe (3) 
It is easily shown that 
dj _ __evn(d) V(bo) 
ate (5), ~ Wnek Tb, exp | - RT ] Gensath, eel (4) 


The ratio |dU/dj|y_) (=R, say) is sometimes referred to as the zero resistance of © 


the rectifier, and its variation with temperature should satisfy (4). Note that — 
with the present model of fixed impurity centres bj, and hence V(b), are 


independent of temperature. 


To discuss the temperature dependence of Vp, take the zero of the potential / 
energy of the current carriers at the top of the full band for p-type semiconductors, — 
and at the bottom of the conduction band for n-type semiconductors, assuming — 


the energy as positive for the states in the almost full band for p-type 
semiconductors, and in the almost empty band for n-type semiconductors. It 
follows that the energy «* of the Fermi level and the energy «’ of the impurity 


levels will in general be negative in both cases. (This is discussed more fully by — 
Landsberg (1952).) Assuming the presence of only one type of impurity | 


(m, centres per unit volume) and neglecting intrinsic conduction (Mott and 
Gurney 1940, p. 175) 


e* ~he’ —EcRT, where c=log (44/7 9/m,) and my =2(27mkT)?2/h8. ..... (5) 
If ¢ be the work function for the current carriers, the diffusion potential V, is — 
elven Dy Vigwdhetced tide: sick Ly ieee eee (6) 
Inserting (5) and (6) into (4) yields 

R, Uy T; tg B 1 1 ‘ ea? e ot 

es = (Ft) exp| F fh on | where B= 6-44 BR Mee (7) 


The variation of c with temperature has here been taken into consideration, and 
suffixes denote values corresponding to different temperatures. B is approxi- 
mately equal to V(4)) and independent of temperature, since experiments on 


Se Re Qin ees 


a ee eee 


reine Soe 


General Theory of Rectifiers 399 


) p- and n-type silicon and on n-type germanium have shown that the change in 


work function as these substances are heated from 200°c to 750°c is quite 
negligible (Smith 1949). It is therefore likely that the corresponding changes 
between room temperature and 200°c are also small. For rectifiers having an 
arbitrary distribution of impurity centres, but subject to diode theory, the equation 
corresponding to (7) is 


R, Ms\H2 (T,\u4 ivy eal 1 
mw Gm) Cr) e[RG-m) | ‘e 
where m, and mz, are the effective masses of the current carriers at temperatures 
sand J... 
Equations (5) and (6), on which these results are based, deserve further 


comment. ‘The temperature variation of Vy has been inferred from forward 
characteristics, and laws which are approximately the same form as (6), i.e. 


Vo=p—q!, g=12-5[log (m)/2,)—O0-121]10-2 evra. -...... (6’) 
have been obtained. ‘The equation for q is the theoretical expression obtained 
from (6) for room temperatures. ‘Table 1 gives experimental results for g, all 


Table 1 Values of g from Forward Characteristics 


gq Xx 10° (ev/°c) Substance Observer 
: Silicon Seiler 1950 
4-7 Selenium Ewels 1950* 
5-9 Selenium Cooper 1952+ 


* The author is indebted to Dr. H. K. Henisch for drawing attention to this reference. 
+ The author is indebted to Dr. R. Cooper for sending him the curve before publication. 


of which would lead to m,~10!%cm-3 if (6’) were used. ‘This equation can, 
however, give only a rough approximation, for the following reasons. (i) Surface 
states would modify the result (6). But there has not so far been any need to 
introduce this notion into the theory of rectification used here. We merely use 
the experimental fact, mentioned above, regarding the negligible temperature 
dependence of the work functions of silicon and germanium (cf. also Markham 
and Miller 1949). (ii) In the case of selenium, internal barriers may complicate 
the situation. (iii) In the case of silicon and germanium, intrinsic conduction 
should be considered. Its effect is to increase the current carrier concentration, 
so that the Fermi level decreases more rapidly with temperature than it would 
in the absence of intrinsic conduction. Actual computations of the temperature 
dependence of the Fermi level, taking into account intrinsic conduction, show that 
theory can lead to slopes of the order of 10-*ev/°c (Hutner, Rittner and du Pré 
1950, figs. 4 and 8, Seiler 1950, fig. 3). ‘Thus, although the inference of Vp from 
forward characteristics is far from certain, experiment and theory are in broad 
agreement that a law of the form (6’) should hold for a limited range of 
temperatures. 

Experiments on the variation of zero resistances with temperature have 
recently been reported (Billig and Ridout 1951), and been interpreted on the 


assumption Rac exp (RM oe a Ra einer! (9) 


in all cases (diffusion and diode theory). As a result Vj was found to increase 
markedly with temperature for selenium, silicon and germanium rectifiers. 
These results, which do not agree with those discussed above, must be due to 
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the use of (9) instead of (7) and (8) in the interpretation of the experiments. 
It may therefore be of interest to apply these latter equations to the experimental | 
curves. 

The expression (7) has been applied to the results for the selenium rectifier 
using 6 = 0-3 ev, and (8) to the results for the germanium rectifier, using B = 0-4 ev, | 
and 17° c was chosen as the standard temperature 7, in both cases. ‘The curves 
have not been extrapolated to lower temperatures, as only one experimental | 
point is available in that region.* The results are given in figs. 1 and2. ‘They 
are to be regarded as semi-quantitative only, owing to the arbitrary choice of , 
the neglect of intrinsic conduction, and the fact that thermal contraction may 
cause complications in the temperature dependence of the zero resistance. 
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Fig. 1. ‘Temperature dependence of the mobility Fig. 2. Temperature dependence 
of the current carriers in selenium. a, from of the effective mass of the 
the zero resistance of a selenium rectifier, current carriers in a germanium 
b, from Hall effect measurements on micro- rectifier from measurements 
crystalline selenium. of the zero resistance. 


No reliable information on the temperature variation of the effective mass 
of current carriers in barrier layers appears to be available, but fig. 1, curve a, 
may be compared with what is known about the mobility of carriers in selenium. 
Mobility determinations from Hall effect measurements are somewhat uncertain, 
but fig. 1 appears to be in rough agreement with the results obtained by Plessner 
(1950) in measurements on single crystals between the temperatures 300° K and 
400°. He found an increase in the mobility by a factor of 5 in this range. The 
existence cf a minimum in the mobility-temperature curve has also been inferred 
by Plessner (1951) from Hall measurements on micro-crystalline selenium, but 
it occurs at a higher temperature (cf. fig. 1, curve 6). Attention has previously 
(1, II) been drawn to the fact that the absolute value of the mobility obtained 
by Henkels (1950) for micro-crystalline selenium, prepared in air, quenched, 
and crystallized at 110°c, was in good agreement with that obtained from 
rectifier characteristics. Henkels found that for this sample the mobility 
increased by a factor of 7:5 as the temperature was raised from 323° to 373°K. 
No independent results appear to be available for selenium in the low temperature 
range. Such direct measurements would be of interest. 

These mobility curves for selenium differ from the corresponding curves 
derived by Pearson and Bardeen (1949) for silicon. Most of their published 
mobility curves have maxima, and their shapes suggest that lattice scattering, which 
is proportional to 7%? (Wilson 1936, p. 211, Seitz 1948), dominates at higher 
temperatures, while impurity scattering, which is proportional to 7? (Conwell 


* It appears that new and improperly understood effects occur at these lower temperatures. 
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and Weisskopf 1950), dominates at lower temperatures. Nevertheless, mobility 
curves of the selenium type have been found fairly frequently (Shockley 1950, 
p. 288) in high resistivity samples of silicon and germanium. Although this is 
not yet understood quantitatively, and selenium has a structure very different 


from that of silicon and germanium, the similarity between these curves should 
be pointed out. 


§3. TEMPERATURE DEPENDENCE OF THE CURRENT-VOLTAGE 
CHARACTERISTICS 


It has recently been found by Cooper (1952)* that two reverse current—voltage 
characteristics of the same selenium rectifier, but taken at different temperatures, 
may intersect. It is perhaps of interest to point out the reason for this effect. 
; Let suffixes 1 and 2 refer to different temperatures, and let 7, be greater 
than 7,. It then follows from (2) and (3) that j, will be greater than, equal to, or 
smaller than j, for a fixed applied voltage (i.e. for a fixed maximum potential 
energy in the barrier) depending on whether 


v1 (Ts) eee Se. 
le exp | 10) Gare? 


Equality can therefore never occur if 7, greater than 7, always implies v, greater 
than v,. In this case the higher temperature always corresponds to the higher 
current for all applied voltages. A necessary condition for the occurrence of 
the cross-over phenomenon is, therefore, that at least part of the mobility— 
temperature curve shall have a negative slope, as it has in selenium (cf. fig. 1). 
If one has two characteristics such that the one for the higher temperature 
corresponds also to the lower mobility, then the two characteristics will intersect 
at an applied voltage given by the above relation, using the equality sign. 

Suppose cross-over occurs when V(b)=V,. Below the cross-over voltage 
V(b)>V,. the left-hand side becomes greater than unity, and 7; >j,. Above the 
cross-over voltage the reverse is the case. It follows that up to the cross-over 
voltage the higher temperature leads to the greater current, while above the 
cross-over voltage the higher temperature corresponds to the lower current. 
These deductions, which are independent of the distribution of impurity centres 
in the barrier layer, explain Cooper’s experimental results. 


§4. THE GENERAL CURRENT-VOLTAGE CHARACTERISTIC 


It will be shown in this section that (2) can be expressed in a different form 
to yield an alternative statement of the theoretical d.c. characteristics which has 
also received attention. It has been found (I, eqn. (13)) that for both a Mott and 
a Schottky barrier, and for the general class of barriers discussed in II, 


ye ce oe EE a (10) 


V, and r are characteristics of a given barrier. For a Mott barrier r=2 and 
V, =edk?T?/é, where « is the dielectric constant of the bulk semiconductor and 
d the thickness of the barrier. Fora Schottky barrier r =4 and V) =ck*T*/87n,e8, 
where 7, is the volume density of ionized impurity centres in the barrier 
(cf. I and II). 

Consider now all barriers for which a relation of type (10) holds. Let the 
reverse voltage required to depress the maximum potential energy in the barrier 


* The author is indebted to Dr. R. Cooper for letting him see these results before publication. 


402 P. T. Landsberg 


to zero be Uy. The condition for this to happen is, by (3), that Vo 
Hence by (2) and (10) Ug satisfies | 


Vee (Vn = URI oy. so ths eee (10’) 
Using (10’) to eliminate Vy from (2), one finds 


of : s (> =" exp eV Ele eee eee (11) 
D 7 


where we have put 


= Vp—U\I 
= 25 fee) oe I Tl aed beeen Al ne eee 12 
is Vo| (=z) | (i 
2RT\12 ven(d) (V pir Ns 
ay ( eae erate ts 1 82° oT Ra 13 
and J ( a ) 2 C= (13) 


The quantities which occur have the following interpretation. Ug is the 
voltage first regarded by Sachs (1942) as the breakdown voltage of the rectifier, 
but its numerical values are too high for this to be the case (II, p. 481 and 
Landsberg 1951c, fig. 2). V can be shown to be the maximum potential 
energy V(b) of the current carriers in the barrier, since by (2), (3), (10) and (10’) 


ip he 
V(b) =Vp—ykRT=Vy—kRT pop Paez) | Si (Peary (14) 

Equation (11) yields, with r=4, the equation used by Billig (1951 b) for a 
Schottky barrier, and with r=2 it yields the equation used by him for a Mott 
barrier, provided a simple manipulation is performed and Vy is neglected 
compared with U,. All characteristics of this class are therefore fully 
equivalent with those obtainable by using Laplace’s approximation (Landsberg 
1949). It should be noted, however, that (11) is somewhat less general than the 
original result (2), since it involves the additional assumption (10). 


§5. THERMAL INSTABILITY OF A BARRIER HAVING AN 
ARBITRARY IMPURITY CENTRE DISTRIBUTION 
Let T be the operating temperature of a rectifier, 0 the ambient temperature, 
so that AT=T—6 is the excess temperature. Then it has been suggested by 
Billig (1951 a, b) that there exists a maximum power loss per unit area per unit 
excess temperature which can be withstood by the rectifier. This value was 
assumed to be a constant C=P/AT for a rectifier of fixed construction and 
properties, which is independent of the ambient temperature. It will be shown 
here that the class of rectifiers for which this assumption is valid can be treated - 
by a uniform method which covers all distributions of impurity centres. 
If the ratio P/AT has its maximum permissible value, and the temperature 
is increased by an amount 57 at constant applied voltage, then the current will 
increase, and the new power loss per unit area per unit excess temperature is 


P d [TP 
AT en Fea 


The condition for thermal instability is that this new value is greater than the 


original (maximum permissible) value P/AT. The condition for instability is 
therefore P 


AT> (@P/OT)y' erccee Cl) 
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For reverse characteristics of rectifiers (2) holds, so that 
eP= Uy" (dittusiom theory)a° 3 ate ne (16) 


Using (15) and (16), one finds that the condition for instability is quite generally 
re VLR AY Le tek ae 


AT = RT 2 var see eee (17) 
‘This shows that the temperature dependence of the mobility may have an important 
bearing on the result. If this effect is neglected, and a Schottky barrier presupposed, 
(17) gives a result first obtained by Billig (1951b, eqn. (21)). For rectifiers having 
_an arbitrary distribution of impurity centres, but subject to the diode theory, the 
condition for thermal instability is different : 


T V(b) 1 TT dm* 


AT ORT 2 Int aT, ee) 
The temperature variation of Vp has been neglected in these equations. Let 
2s ce for the diffusion theory, 
T Oar 
Beers  UewhereD= 3 gh ee (19) 
aM —_ Om" for the diode the 
— ape gp for the diode theory, 


then conditions (17) and (18) mean that thermal instability occurs if the 
temperature rise is large enough, so that y<V(b)/RT. ‘Thus the Sachs criterion 
V(b) =0 (cf. §4) is modified to V(b) =ykT 


in the present model. Let the values of the parameters be distinguished by a 
suffix B if this condition is fulfilled. It then follows from (14) and (20) that 
Vp WGN neCAT 
et ae) Ds Up= sees es 24 
Pan T yn FE ila) ik 


Billig determined the value of Uy by a graphical method, namely, by finding 
the turning points of lines of constant C which cut across theoretical current— 
voltage curves for different temperatures. For both diode and diffusion theory, 
and for all impurity centre distributions, the value of Ug so determined is the 
same as that given by (17) and (18) (cf. Appendix). It will now be shown that the 
theory can be carried through analytically. 

It is helpful to make the assumption (10) so that when |U|>Vj, 
V)=(2cb/e)"U(RT)". As Vo is independent of applied voltage, and 6 is 
independent of temperature, it follows that the potential energy V, defined by 

Va (Gy Seb) Ce it, AS oe oh. (22) 


is independent of voltage and temperature. Its values are given in table 2. 


Table 2. Values of V, (Notation as in IT) 
e 


y=2 Mott barrier Vi= 


v0) fe 
Pa 


Schottky barrier Barrier introduced in II 
V3=87re'n,/<8 V3 =16e no/e® 


6 

General class y= 47e 
. Lit 3 
of barriers € 


404 P. T. Landsberg | 


Assumption (10) leads with the notation (22) to a second equation for Ug: 
Uz w) kT3 s—1 - 
RT, = ( V, ) YVR- Pry te ere as (23) 


(22) and (23) are the two equations for the temperature 7, and the applied — 
voltage U, at which instability sets in. The temperature must be determined 
first, and is found in both diffusion and diode theory to be given by 


(yak Tayi Vale NT ae yee ee (24) 


Introduce the notation 
B=Vp )/k6, 7=AT/0, 
- 5)" aNRY Ga) for the diffusion theory, 
Cé 
n(d)(k0)?" 
Substituting in (21) and (24) from the general theory of barrier layers (cf. II), 


and using (25), it is found that the maximum permissible temperature rise is 
given by 


f B r+3/2 
mu teye| rt 94D | exp (—7 1—4+D) 


(2nm*) 42 


Vay pe eee 
ko or e diode cory. 


Oe for the diffusion theory, ...... (26 a) 
con NE ie a oe E =. as +D | exp (ieee 
for the diode theory. ...... (26d) 
The breakdown voltage is given by 
poe {( ) Fe See D | \ in both theories (26 c) 
(kay ES i fg eee 


For a theoretical calculation it is necessary that A, B, D, RO and r be known. 
It would be useful to supplement eqns. (26) by graphical information, but this 
is postponed since the basic physical ideas involved are still under discussion 
(Cooper 1952). 

Only one application of the equations (26) will be given. The three functions 
on the right-hand sides increase with 7 unless D is negative, when they may 
attain a maximum and subsequently decrease. If the thermal dissipation 
constant is decreased, A decreases, hence + decreases, and therefore the 
breakdown voltage must decrease, provided D is not negative. Cooper, however, 
observed an increase of the breakdown voltage when the thermal dissipation 
constant was decreased. If a theory based on thermal instability alone can 
account for these observations at all, then only if D can be shown to have a negative’ 
value. ‘This shows, incidentally, that in any test of the theory it is important 
to take into account the temperature variation of the mobility if diffusion theory 
is used (or of the effective mass if the diode theory is used). 


§6. CURVES OF CONSTANT POWER DISSIPATION 
The power dissipation per unit area P and the maximum permissible 
temperature rise AT depend on several parameters: C, 0, «, Vp, etc. Suppose 
any one of these, p* say, is allowed to vary while all the others remain fixed, 
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and let suffixes 1 and 2 refer to the conditions when p* has the values p,* and 


) § p.* respectively. ‘Then, according to the assumption stated at the beginning 


of §5 (to be referred to as A), the voltages U,/|e] and U,/|e] at which thermal 
instability sets in are related by P(U,, p,*)/AT, = P(U,, p,*)/AT,. If, however, 
the assumption of §5 is not valid, and it is more correct to assume that a certain 
rectifying junction can withstand only a certain maximum power loss P per unit 
area (this assumption is to be referred to as B), then U, and U, are related by 
P(U,, ~,*)=P(U,, po* =P. This condition yields curves of the breakdown 
voltage against p* on the basis of assumption B (which must be clearly 
distinguished from assumption A). 

Curves of this type have first been discussed for a Schottky barrier by Billig 
(1951 a, figs. 1 and 2, or 1951b, figs. 7 and 8). At an ambient temperature of 
17° c the value of P for the onset of thermal instability was taken as 0-92 watts/cm?. 
The variable parameter p* was the conductivity and the diffusion potential.+ 
The constants used in these curves imply, when evaluated, a dielectric constant 
of about 19 and an impurity centre density » of about 2 to3 x 10!®cm-%. As these 
values are not always likely, additional curves have been calculated (figs. 3 and 4), 
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Fig. 3. Applied voltage against conductivity Fig. 4. Applied voltage against (a) diffusion 
for constant power dissipation per unit potential and (d) dielectric constant for 
area in a rectifier having a Schottky constant power dissipation per unit area 
barrier and subject to diffusion theory. ina rectifier having a Schottky barrier 


and subject to diffusion theory. 


assuming a Schottky barrier and using the values (cf. I, table 2) «=6:3 (Se), 
n=1-3 x10'8cm-3, RT =1/40ev, Vyp=0-5ev, P=0-92 watt/cm*?. The last three 
values were also used by Billig. In fig. 3 a value of 7-1 x 108 volt/cm? sec was 
used for the mobility (cf. I, table 2). 

Another curve of constant power dissipation, based on the same constants, 
is added in fig. 4. It shows the important effect which the dielectric constant 
has on the power dissipation. 

These curves are to be regarded merely as curves of the variable parameter 
against the applied voltage for constant power dissipation per unit area. [or 
rectifiers for which the value of P which is required for the onset of thermal 
instability is 0-92 watt/cm? the applied voltage plotted must, however, be very 
near the breakdown voltage. 


§7. THEORY OF CURRENT CREEP 
Study of the time dependence of the current-voltage characteristic of 
selenium rectifiers (Henisch and Ewels 1950) shows that negative creep—the 
decrease of current with lapse of time for fixed external conditions—is met, with 


+ There is an error in the equation on which this curve was based (Billig 1951 a, eqn. (6), 1951 b, 
eqn. (24)); the ratio of the diffusion potentials should be replaced by unity ; this has a small effect 
on the curve given. 
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quite frequently. Positive creep can occur in selenium rectifiers for very sigan 
time intervals after the application of a voltage, and after special treatment of the 


rectifier a short time before the external potential is applied. Rather different 


phenomena occur for other rectifiers. Copper oxide rectifiers, for instance, 


show mainly positive creep. It is therefore desirable to go as far as possible with — 


the general theory, and thereafter attempt special theories based on specific 
models. 


(i) General Theory 
Since creep occurs mainly for reverse potentials, the current density is given by 
juert?, 
where y is given by (2) and 
avs oe for the diffusion theory (cf. II, eqn. (10)), 
jo for the diode theory (cf. II, eqn. (13)). 
At constant temperature the current density (27) is a function of the time and of 


the applied voltage because the space charge distribution, and hence y, will 
in general depend on these quantities. ‘The greatest part of the time and voltage 


dependence of the characteristic is therefore determined by exp y in both diode” 


and diffusion theory. Since V’(d) is usually negligible, and the plane x=d is 
very near the rectifying electrode, we have approximately from (1), (2) and (27) 


log J +-G4/O5, oa ee (28) 
ae and O=|e| iE n(x) dx. 


Q is the total (numerical) space eb per unit cross section of the barrier; its 
sign will in general be opposite to that of the predominant current carriers. 
As previously suggested (Landsberg 1951a), therefore, the problem of predicting 
current creep 1s equivalent to the problem of determining the time dependence of the 
total space charge Q in the barrier. ‘ Positive’ creep implies that O increases with 
time, ‘negative’ creep that it decreases with time. ‘To determine O(¢), however, 
one must know how the thickness d of the barrier varies with time. Some 
information may be deduced from II, eqn. (3) (cf. also Landsberg 1951d, 
Appendix B), provided the space charge distribution function, including its 
dependence on time, is already known. 

One can go one step further in a general theory, since (Rose and Schmidt 1947) 


where Ge 


{ Fn RAE ERS fice (29) 


must remain constant for fixed external conditions. It follows that if n(x, 2) 


| 
. 
| 
| 
| 


decreases with lapse of time at all points in the barrier, the barrier becomes — 


thicker, while if (x, t) increases, d decreases. In particular, if m depends on 
« and ¢ through two separate functions, n(x, t)=/(x)g(t) say, the function g can 
be eliminated, since it follows that 


fx) g fo Le) 
n(x, t) = S ————._ and 
( b [rate is Q= le |S [xf (x) ae me oe 
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Both quantities depend on time through the barrier thickness. If of (ey xe, 

s+2]e|S 

= a i een 1). ee ee (31) 

This shows the type of relation to be expected between Q and d. Thus, in 

general, positive creep implies decreasing barrier thickness, negative creep 
increasing barrier thickness. 

If no new potential barriers are formed, and if parameters such as Vp and T 
remain unaltered, then the above results are largely independent of the physical 
or chemical processes which may be responsible for current creep or forming. 
They show that these processes are important because they lead to a modification 
of the total space charge Q(t) in the barrier. That is why the above considerations 
may properly be regarded as part of what we have defined (in IT) as the general 
theory of rectification. 


(ii) Special Theories 

No quantitative special theories are so far available. But they will have to 
cover at least two physically important cases. First, the case in which a new 
compound is formed at the contact. The theories will here have to deal with the 
passage of current carriers (electrons or ions) of one type from the bulk 
semiconductor, while those of the opposite charge migrate from the metal, or, 
possibly, from the material trapped between metal and semiconductor. ‘The 
carriers presumably combined near the contact to form a new phase of growing 
thickness (e.g. a selenide). The time rate of growth of this phase is clearly 
important, and governs the time-dependence of the current. There may also 
be a limiting thickness beyond which the new substance cannot grow. Processes 
of this kind may be expected to offer analogies with the growth of oxide films 
on metals, and the special theories will therefore have to be based on what we know 
already about tarnishing reactions. If the new phase is relatively free from 
space charge, there would be a tendency for the barrier thickness to increase 
(by (29)). The total space charge Q in the barrier would drop due to the formation 
of the new compound, and increase due to an increasing barrier thickness. The 
general theory shows that the object of any special theory must be to estimate the 
importance of these two opposing trends, and hence to obtain the time 
dependence of QO. , 

The second class of physically important cases arises if the changes in the 
barrier can be regarded as due to forced diffusion (1.e. diffusion in the presence 
of a field). A tentative suggestion regarding the mechanism of current creep 
according to theories of this type will be offered. In it the assumption is made 
that the diffusion constant incorporates the effect of the electric field. ‘This is 
an over-simplification, but it enables one to obtain a first insight into the problems 
involved. 

Suppose, with Cooper (1950), that positive ions can enter the barrier in the 
case of a p-type semiconductor. If they do not penetrate very far, the effect of 
their presence on the barrier thickness may by (29) be neglected. Since the total 
space charge in the barrier is negative, these ions will, however, cause a decrease 
in the total space charge (in the semiconductor), even if they are eventually 

neutralized by electrons, provided only that these electrons come from the 
semiconductor. Ordinary diffusion theory suggests that a law of the type 
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OQ «(D,t)-¥? may be valid, D, being an appropriate diffusion coefficient, so that 
(28) becomes logj = AOL 0” eee (32) t 
Equation (32) is in agreement with creep experiments on selenium rectifiers,> 


as clear from fig. 5. The quantity B was found to have the value 0:44 min” for — 
the upper and 0-49 min"* for the lower curve. 
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Fig. 5. Current-time curves for selenium rectifiers. 
a. Specimen B2e1, U=16-6 ev, T=22° c. 
b. Specimen B2b2, U=12-6 ev, T=19° c. 


APPENDIX 


Proof of the Equivalence of the Graphical and the Analytical Method of Determining 
the Onset of Thermal Instability 

Consider a barrier having an arbitrary distribution of impurity centres. 
Diffusion and diode theory can be treated together, since the temperature 
dependence of the current is in both cases given by (cf. II, eqns. (11) and (13)) 

j=aT? eT, where —RB=Vp—kTy=V(b). ...... (A 1) 

a is temperature dependent since it involves the mobility if the diffusion theory 
applies, and the effective mass if the diode theory applies: 


const. x v (diffusion theory), 
| const. x 1/m* (diode theory). 


Consider a curve of constant power dissipation C[=jU/e(T—6)] per unit area 
per unit excess temperature, which cuts across current-voltage characteristics 
for various temperatures. For such a curve 


oU eC 27? 
Gy = ce 2B+(2T adalat) | 


since temperature and current are related by (A1). 

Suppose (dU/d/)¢ is negative, i.e. the value of C can be maintained after a 
small rise in current only if U is correspondingly decreased. Then for phenomena 
which occur at constant applied voltage a small rise in current must lead to a 
corresponding rise in C, If C has already its maximum permissible value, 
however, a slight increase in current must lead to thermal instability. Similarly, 
if (0U/d])¢ is positive, conditions are thermally stable. One finds after a short 


transformation that 
oU te ih VN bc 
ect < 0 Peale 2 TH . 
Fj Jo 
} The author is indebted to Drs. H. K. Henisch and J. Ewels for putting the details of their 
experimental results at the author’s disposal. 
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This leads to the conditions (17) and (18) of the text, thus proving the 


equivalence of the graphical and the analytical method for all impurity centre 
distributions. 
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ABSTRACT. The paper describes an experimental investigation of breakdown in 
selenium rectifiers. Cooling from about 110°c to —74°c caused the breakdown voltage to 
decrease continuously, and at the same time anomalous behaviour was observed in the 
current-voltage characteristics. When high voltages were applied, the current decreased 
as the temperature increased, but at low voltages rise in temperature caused the current to 
increase. 

The effect of changing the thermal dissipation constant was similar to that caused by a 
change in ambient temperature. Reduction in the dissipation constant, effected by placing 
the rectifier in a vacuum, was observed to increase the breakdown voltage, but the leakage 
‘current at low voltages also increased. 

None of the available theories appears to account satisfactorily for the observations. 


§1. INTRODUCTION 
LL contact rectifiers break down when excessive reverse voltages are applied. 
The phenomenon is illustrated by the current-voltage (/, V) characteristics 
shown in fig. 1. As the reverse voltage is increased, the differential 
resistance 0V/d/ decreases and finally becomes approximately equal to the resistance 
of the rectifier in the conducting direction. The barrier layer has then become 
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ineffective, and by extrapolating the reverse characteristic in the manner shown in 
fig. 1 the breakdown voltage (U,) may be determined. Breakdown in contact 
rectifiers must not be confused with insulation breakdown; it is not accompanied | 
by disruption of the material between the two electrodes, and it is reversible. 

Although rectifier breakdown is a phenomenon of some practical interest, there 
appears to be a dearth of detailed experimental facts about it. However, theories 
have been proposed by Sachs (1942), Schottky (1942) and Billig (1951), and their 
predictions are compared in this paper with the results of some experiments 
recently carried out by the author. 


it 


Width of loop w : 100v 
increases with V 
frequency 
-74°C 
Fig. 1. Dynamic characteristic illustrating Fig. 2. Illustrating effect of ambient 
breakdown of a selenium rectifier. temperature on breakdown. 


§2. EXPERIMENTAL 

The experiments were carried out to determine how the breakdown voltage 
varied with ambient temperature, counter-electrode metal, and the thermal 
dissipation constant of the rectifier. Both commerical and fully formed laboratory 
prepared rectifiers were used. A vacuum evaporation method similar to that 
recently described by Billig (1951 b) was adopted for the preparation of the latter 
specimens, and spectroscopically pure materials were used. 

Since Joule heating and current creep rendered it difficult to obtain reliable 
static //V characteristics over an extended range of reverse voltages, breakdown 
was investigated by displaying dynamic characteristics on a d.c. oscilloscope. The 
reverse characteristics were of the form shown diagrammatically in fig. 1. The 
width of the two loops depended on the frequency of a.c. supplied to the measuring 
circuit. At 50 c/s and less the loops were not pronounced, but with increase in 
frequency the major loop increased considerably in width. Attention has 
previously been drawn by Billig (1951 b) to the existence of the major loop, but he 
apparently did not study the effect of frequency on it. In the present paper the 
values of breakdown voltages were obtained at 50c/s from the portion of the 
characteristic corresponding to increasing values of reverse voltage. 

‘The ambient temperature was varied from about 30°c to — 74°c in experiments 
with the laboratory prepared rectifiers; in the case of the commercial rectifiers it 
was varied from 120°c to —74°c. At room temperature and above the rectifiers 
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were immersed in transformer oil. Pentane was used below 20°c. The breakdown 
voltage of both kinds of rectifier was found to decrease continuously with fall in 
temperature. ‘The effect observed with a commercial rectifier is illustrated by the 
oscillograms shown in fig. 2. At 110°c the breakdown voltage was about 70 volts, 
at 25°c it was 50 volts, but at —74°c it was only 25 volts. The effect was repro- 
ducible, and immersion in oil at 110°c did not destroy the rectifier, despite the fact 
that the high temperature caused the counter-electrode to melt. 

There is also apparent in fig. 2 another phenomenon of considerable interest. 
Comparison of the characteristics shows that whilst a rise in ambient temperature 
had caused the leakage current corresponding to small reverse voltages to increase, 
the reverse effect had occurred at voltages on the knees of the characteristics and 
beyond. With increase in reverse voltage the temperature coefficient (0/07 at 
constant voltage) had changed from positive to negative. This phenomenon 
occurred with both commercial and laboratory prepared rectifiers. Its existence 
was confirmed by applying various continuous voltages to the rectifiers and at each 
voltage, when the rate of current creep had become small, observing the change in 
current caused by a rise in temperature. 

The characteristics shown in fig. 2 are consistent with the two curves published 
by Komar (1947) illustrating the effect of cooling a commercial selenium rectifier 
from room temperature to about —50°c, but they are not consistent with the 
characteristics given by Williams and Thompson (1941). The latter show that 
heating a rectifier from 20°c to 80°c caused the leakage current to increase at all 
voltages. Possible differences in the manufacturing processes of the rectifiers, 
about which little is known, may account for the difference in behaviour. A negative 
temperature coefficient has been observed previously in laboratory prepared 
rectifiers by Schmidt (1941), who obtained static reverse (J, V) characteristics up 
to about 1-0 volt at temperatures between 22°c and —92°c, and also by Henisch 
and Ewels (1950), who found that cooling a rectifier from room temperature to 
about —50°c with 10 volts d.c. applied caused the leakage current to increase. 

The effects just described were found to be independent of the metal used as 
counter-electrode. But at room temperature the breakdown voltage differed 
considerably for different metals. Experiments were made using rectifiers with 
counter-electrodes of magnesium, aluminium, cadmium, bismuth, zinc and 
antimony. Variation was observed in the properties of apparently similar rectifiers, 
but the effect of counter-electrode metal on the breakdown voltage at room 
temperature is illustrated by the following average values : magnesium 130v, 
aluminium 100 v, cadmium 60 v, bismuth 40 v, zinc and antimony about 20 v. 

Ina recent paper (Billig 1951 b) it has been stated that the breakdown voltage of 
selenium rectifiers is increased by increasing the thermal dissipation constant. 
This statement was not borne out by the following experiment. A commercial 
rectifier was suspended in a bell jar, and by means of a suitably placed valve a 
stream of air could be drawn over the rectifier; alternatively the bell jar could be 
evacuated to a pressure of about 10->mm Hg, thus changing considerably the 
thermal dissipation constant. 

In fig. 3, curve a is the dynamic characteristic observed when the air stream 
passed over the rectifier. ‘The temperature of its counter-electrode, measured by a 
thermocouple, was 28°c. Curve 6 is the dynamic characteristic observed after 
the air pressure in the bell jar had been reduced to about 10-° mm of Hg. During 
the time taken to reach this pressure the voltage applied to the measuring circuit 
was kept constant, but the temperature of the rectifier rose above the melting 
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point of the counter-electrode, approximately 120°c, before thermal equilibrium 
was restored. ‘The effect brought about by the considerable reduction in the 
thermal dissipation constant is clearly similar to that caused by increasing the 
ambient temperature from about 25°c to 100°c, and illustrated in fig. 2. The 
leakage current increased at low voltages, but breakdown occurred at a higher 
voltage. 


@& Commercial Rectifier 


Vq (volts) 


4 Antimony 
10->mm H Counterelectrode 


‘Temperature 120°C 


Atmospheric Pressure 200 250 300 350 b 


Temperature 28°C Temperature (°k) 
Fig. 3. Effect of change in thermal Fig. 4. Effect of temperature on Vg. Result obtained 
dissipation constant on breakdown. from forward (J, V) characteristics. 


§3. DISCUSSION 
3.1. Image Force Theory 

Since the electric intensity in a barrier layer is estimated at about 10° or 10° v/cm, 
the height of the potential barrier may be influenced by image force, and Sachs 
(1942) has suggested that breakdown will occur when the applied voltage is sufficient 
to depress the height of the barrier to about the same magnitude as the average 
thermal energy of current carriers in the semiconductor. When this criterion 
is applied to the Schottky type of barrier layer, which is generally assumed to 
exist in selenium rectifiers, the following equation results : 

Ke Vt 
O68 NT ee ee (1) 
where U,=breakdown voltage, K = permittivity of the semiconductor, e= charge 
onan ionized impurity centre, N = density of charged impurity centres, V, = height - 
of potential barrier negiecting image force. 

Landsberg (1951) has previously pointed out that when reasonable values for 
the parameters are inserted the value obtained for U, is several orders of magnitude 
too great, but eqn. (1) is based on a very simple model of the barrier layer that might 
be expected to yield high values of U,. ‘The theory cannot be rejected for the 
above reason alone. However, the theory appears to be contradicted by the 
experiments since it apparently predicts U, independent of temperature. But 
since both V4 and N may vary with temperature, the effect of possible variations 
in these quantities must be considered. 

Equation (1) indicates that U, is most sensitive to variations in V,, and Billig 
and Ridout (1951) have deduced from the variation with temperature of the zero 
resistance (0V/01),_,, that V rapidly decreases with fall in temperature between 
20°c and —195°c. However, this result has been questioned by Landsberg 
(1952) who has shown that these measurements can be interpreted quite plausibly 
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by assuming constant V, and temperature dependent current carrier mobility. 
) Clearly the conclusion of Billig and Ridout cannot be accepted unlessit is confirmed 
by an alternative method of measurement. 

The variation of V with temperature can also be obtained from Schottky 
plots of (1/capacitance)? against bias voltage, and also by extrapolating the linear 
portions of the forward (/, V) characteristics back to the voltage axis. The first 
method is, however, not a good one. It involves extrapolating curves which, in 
experiments made by the author, were generally not linear ; also the capacitance of 
selenium rectifiers is in many cases (Cooper 1949, Joos 1948) a function of frequency, 
and until a satisfactory explanation of the effect is found some doubt must exist 
about values of V, and N obtained by this method. However, the author is 
unaware of any capacitance measurements which indicate that V., decreases with 
fallin temperature. ‘The second method has the advantages that the observations 
require for their interpretation no formula based on a particular model of the 
barrierlayer. Figure 4 illustrates results obtained by this method on a commercial 
rectifier (curve a) and on a laboratory prepared rectifier with a counter-electrode of 
antimony (curve 6). In both cases Vq increased linearly with fall of temperature 
from 20°c to about —74°c. This result is qualitatively in agreement with the 
theory of Fowler (1936), but numerical agreement is not good. The contact 
potential difference between a metal and an exhaustion semiconductor like 
selenium with about 10!’impurity centres per cm® should, acccording to the theory, 
increase linearly only by about 0-05 ev on reducing the temperature between the 
above limits. ‘The reason for this difference is not known, but errors caused by 
effects at grain boundaries may exist in values of V, obtained from forward 
characteristics. It is, however, unlikely that they alone would account for the 
difference in results obtained by the two methods. 

For the above reasons it cannot be assumed that V, decreases with fall in 
temperature, and an appreciable increase in N must be assumed to account for the 
observed decrease in breakdown voltage. But capacitance measurements made by 
Schmidt (1941) provide no evidence of this, nor have similar measurements made 
by the author. However, Plessner (1951) has given reasons for believing that the 
density of impurity centres in micro-crystalline selentum may be increased by 
reducing its temperature below 20°c. But it is not certain that the material of 
the barrier layer possesses the same properties as the bulk Se: there is some 
evidence for believing a chemical compound is formed at the interface (Williams 
and Thompson 1941, Henisch and Ewels 1950). 

It does not appear that N increases with fall in temperature, but this 
possibility cannot be definitely ruled out because of the sparse and unsatisfactory 
available evidence. 

The experiments described in §2 have shown that at room temperature the 
breakdown voltage depends on the counter-electrode metal, and this result is 
qualitatively in accordance with the image force theory of breakdown. 


3.2. Schottky’s Pass Effect 


Capacitance measurements on selenium rectifiers indicate that the barrier 
layer thickness is only equal to a few times the average separation between charged 
impurity centres. This led Schottky (1942) to point out that we cannot assume a 
continous variation in potential through the barrier layer. We must assume that 
it consists of a range of hills separating the neighbouring charged impurity centres. 
According to Schottky, enhanced current at low voltages and breakdown at higher 
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voltages may result from the existence of localities, called ‘passes’, that 
fortuitously contain a greater density of charged impurity centres than the average 
for the barrier layer. In effect a ‘ pass’ provides a path of low resistence through 
the barrier layer. ‘The magnitude of the resistance may not greatly exceed the 
spreading resistance of the semiconductor in series with the pass. 

The theory does not allow calculation of the breakdown voltage, but a decrease 
in this quantity with increase in N may be expected. In this respect the theory is. 
similar to the image force theory, and the possibility of NV increasing with fall in 
temperature has been discussed. 


3.3. Thermal Instability 


Breakdown in rectifiers may occur (Billig 1951 a, b), provided the temperature 
coefficient is positive, when the heat generated per sec per °C above ambient 
exceeds that dissipated per sec per °c above ambient by the cooling arrangement. 
But the rectifiers used for the experiments described in this paper possessed a 
negative temperature coefficient at voltages on and beyond the knee of the (/, V) 
characteristic, and for this reason breakdown cannot be attributed to thermal 
instability. For the same reason the loop in the breakdown portion of the 
dynamic (J, V) characteristics cannot be accounted for in the manner suggested 
by Billig (1951 b). 

§4. CONCLUSIONS 

‘The experiments have shown that breakdown may occur in selenium rectifiers 
at higher voltages as the ambient temperature is raised, and as the thermal 
dissipation constant is lowered. At low voltages the temperature coefficient is 
positive, but near to and beyond the knee of the (J, V) characteristic it is negative. 

The results do not support the thermal instability theory of breakdown. The 
other theories cannot be ruled out because not enough is known about the variation 
with temperature of the density of charged impurity centres in the barrier layers. 
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“ABSTRACT. Experimental studies have been made of the distribution in size of output 
pulses from a scintillation counter detecting monoenergetic alpha-particles. Separate 
investigations have been made of contributions to the pulse size spread by the photomulti- 
plier and by the scintillating crystal. It is shown that the main effect is due to statistical 
fluctuations in the number of photoelectrons produced at the photomultiplier cathode. 
The effect of the multiplying stages on these fluctuations is investigated and also the effect 
of flaws in the detecting crystal. The frequency distribution of the pulses is usually found. 
to follow the gaussian form. 


§1. GENERAL INTRODUCTION 


HEN ascintillation counter, comprising photomultiplier and luminescent 

crystal, is used to detect monoenergetic particles fluctuations occur in 

the magnitude of the output pulses. Several workers have made 
experimental and theoretical studies of such fluctuations (Engstrom 1947, 
Morton and Mitchell 1948, Morton 1949, Hoyt 1949). In order to determine 
the various sources of fluctuations, a procedure not adequately covered by the 
above workers, we have made separate studies of the contributions from the 
various stages of the photomultiplier and from the phosphor crystal under the 
conditions obtaining in scintillation counters. ‘The following factors are likely 
to contribute to the pulse spread: (i) variation in intensity of successive light 
pulses arriving at the photomultiplier cathode due to local variations in light 
production in the crystal and to variations in light transmission from crystal to 
cathode, the latter arising from crystal flaws, reflection and absorption; 
(ii) fluctuations in the number of electrons passing through the photomultiplier 
due to low photoelectric response of the cathode and non-uniform response 
over its surface, to low electron collection efficiency in the multiplying stages 
and to variation in response over the dynode surfaces. 

In order to investigate effects due to the photomultiplier the phosphor crystal 
was replaced by a neon lamp connected in a simple saw-tooth oscillator circuit. 
The light pulses from this source are found to be of uniform intensity and have 
a controllable repetition frequency. By reducing the light intensity with neutral 
filters, calibrated by use of a photometer with similar neon lamp source, the 
light pulse intensity at the photomultiplier cathode can be varied over a wide 
range. After investigation of photomultiplier characteristics, measurements 
with the phosphor crystal as light source enable the contribution of the latter 
to pulse variations to be determined. (Since the neon lamp light pulses are of 
relatively long duration, about 10sec, compared with scintillation pulses, the 
time constant of the multiplier output circuit is made large enough to integrate 
the pulses and so provide a true measure of total light intensity per pulse.) 

2 Bas 


416 G. F. F. Garlick and G. T. Wright 


§2. THEORETICAL INTRODUCTION 


In other investigations, referred to above, a mean quantum efficiency for the 
whole of the photomultiplier cathode has been assumed and also a mean gain 
for each multiplying stage. In a detailed analysis, not included here, we have — 
found that such average values are justified if the average photon distribution 
over the cathode surface remains the same for successive scintillations and if the 
average distribution of incident electron density on any dynode surface is the 
same for successive electron pulses. There is an important difference between 
our studies and previous investigations. Statistical considerations of other 
workers are developed by determining the output pulse magnitude due to single 
electrons leaving the photocathode, whereas in our studies the practical conditions 
are assumed. We therefore consider output pulse fluctuations due to any given 
number of photoelectrons leaving the cathode. A further difference is that 
Hoyt (1949) considers separately fluctuations due to the multiplying stages but 
considers collectively effects due to the photocathode, light collection system 
and phosphor crystal. It must also be pointed out that the considerable 
discrepancy found by Morton and Mitchell (1948), and Morton (1949), between 
experimental and theoretical results is due partially to a mistake in their theoretical 
derivations. To correct for this discrepancy they assume interstage focusing | 
losses for electrons. However, in their calculations they use the secondary 
electron emission ratio of the dynodes in an expression where the mean gain 
per stage is what is required. As a result by later correcting for focusing losses ~ 
they include the effect of such losses twice. 

Following analyses similar to those of the above workers, we obtain an 
expression for the variance of the charge arriving at the collector electrode of the 
photomultiplier per scintillation as follows: 

GO: 
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where barred symbols indicate mean values, e is the electronic charge, X the 
number of photons produced by absorption of one monoenergetic particle in 
the phosphor crystal, p is the probability that a photon incident on the cathode 
produces an electron which enters the multiplying stages, f is the fraction of the 
photons produced which reach the photocathode, s is the mean gain per 
multiplying stage and ¢ the number of stages. By assuming s+! 1 and taking 
M as the mean gain of the multiplying stages we have 

OG — O; yo  y2 2 VY Ss 

as = (X2—X")p?f2 M? + pfX M? (5 -?f). ssa rartt (2) 
The first term on the right-hand side represents the contribution from the 
phosphor crystal. ‘The second term is similar to that of Morton and Mitchell 
for the effect of the photomultiplier but it now includes the product pf in 
parenthesis which takes account of the efficiency with which photons originating 
in the crystal produce electrons which enter the multiplier. For the two types 
of photomultiplier examined (EMI types 4588 and 5311) and their accompanying 
crystal systems pf is about 0-01 and can be neglected compared with the factor 
s/(s—1). Its small size is due to low light collection efficiency in the case of the 
arrangement using the 4588 type tube and to the loss of light by transmission 
through, and reflection from the cathode of the 5311 type multiplier. 
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When pulses of electrons leave the photomultiplier cathode the resulting 
pulse distributions are different from those due to single electrons leaving the 
cathode (cf. Engstrom 1947). In our experiments we find the distributions to 
be of gaussian form with a given half width 7. The-use of the usual gaussian 
formulae together with eqn. (2) gives the following expression for 7: 


= B+ (a), iat (3) 


where K=5-:54 eM/c, c being the capacity to earth of the collector electrode, 
V is the mean pulse voltage and A(£) the contribution of the phosphor crystal. 


§3. EXPERIMENTAL STUDIES 
(i) Experimental Methods 


‘Integral bias curves’ were obtained by feeding the output pulses from 
the photomultiplier into a conventional discriminator and scaler unit. 'To obtain 
the derivative of these curves directly a single channel pulse analyser of unique 
design was used. ‘This apparatus, to be described elsewhere, presents the 
required frequency distribution as a trace on a cathode-ray tube screen which 
can be photographed. 

When using a luminescent crystal and alpha-particle source to provide 
scintillations the efficiency of light collection was improved by cementing the 
given small crystal to the glass envelope of the photomultiplier. In the 
experiments to be described an EMI type 4588 multiplier was used with 
opaque antimony-—caesium photocathode (side window entry for light) and 
nine stages of multiplication. However, similar results have been obtained for 
an EMI type 5311 multiplier. All voltage supplies were very carefully stabilized 
and very slow d.c. drift of the supply was also compensated. The latter is a 
serious disturbing factor when making observations over several minutes with 
a conventional scaling unit. 


(ii) The Form of the Pulse Distribution Curve 


If the pulse distribution is of gaussian form, then if V is the mean pulse height 
and V any given pulse height the half width of the distribution 7 is given by 
n=(5:540)!?/V = —2-36/mV, where o is the variance of the distribution and 
m=y/(V—V), y being a function which is obtained from tables (Jahnke and 
Emde 1933), using the experimental values of counts per second for a given 
pulse height V. In fig. 1 we give the relation between y and V obtained in this 
manner from ‘integral bias curves’ by using various luminescent crystals and a 
homogeneous alpha-particle source (5 Mev) and also by using the neon lamp source. 
It can be seen that the experimental values confirm the expected linear relation 
between y and V. The broken horizontal lines in fig. 1 indicate the values of 
y for 98% (A) and 2% (B) of the total counts per second. Experimental 
observations outside these limits are not useful for obtaining values of y, as the 
relative statistical errors become large and are much magnified in the final value 
of y. Curves 1, 2, and 3 were obtained with thin flat crystals of scheelite, 
anthracene and naphthalene containing stilbene. Curves 4 and 5 were for 
crystals with many internal flaws and surface irregularities, while curve 6 was 


for the neon lamp pulses. 
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(iii) Fluctuations in the Photomultipler 


By using the neon lamp source we eliminate pulse spread due to the photon 
source and so we can investigate the dependence of pulse distribution on | 
multiplier characteristics. q 

(a) Effect of multiplier gain. Using a low intensity steady light source the | 
gain of the multiplier was investigated under ‘static’ conditions by comparing 
photocathode and collector electrode currents. It is found that the gain M is 
given in terms of the overall applied voltage v (equal volts per stage) by 


M = of ago." on Ta ane ee (4) 


For the multiplier used vp =98-2 volts and y=5-35 respectively. By replacing : 
the steady source with an anthracene crystal and alpha source the gain was 
investigated under pulsed conditions. Equation (4) was found to hold provided 


VIN ARB. UNITS = ———> 


Fig. 1. Variation of the function y of eqn. (8) with scintillation pulse magnitude. 


1. Scheelite. 2. Anthracene. 3. Stilbene in naphthalene. 4. Copper activated zine 
sulphide. 5, Thallium activated potassium iodide. 6. Neon lamp pulses. 


that the interstage capacity was made large enough to avoid voltage changes 
between dynodes during the transit of a pulse of charge. A value of y=5-34 
was obtained over a range extending to 1500 volts (150v/stage). At 1500 volts 
the gain was 2:2x10°. The dependence of the half width y of the pulse 
distribution on gain was investigated and was found to be small, as would be 
expected, since from eqn. (3) 7=%o{s/(s—1)}?._ With the above multiplier 
constants the measured value of 7 changes from 0-19 at 700 volts (overall) to 
about 0:17 at 1200 volts for the given mean pulse magnitude. From 
measurement of the effective collector electrode capacity and the overall gain 
it was deduced that an average of 360 photoelectrons were produced per 
scintillation. If we assume that all photoelectrons enter the multiplying stages 
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and also apply the correction factor {s/(s— 1)}1? then the half width expected from 
statistical fluctuations in this number is 0-135. Since the fraction of photo- 
electrons entering the multiplying stages and the correction factor {s/(s—1)}? 
are inter-dependent quantities, the true value of each must be obtained by 
successive approximations. However, a change from 1-00 to 0-60 in the fraction 
of photoelectrons collected alters the quantity {s/(s—1)}4? by less than 1% and 
so only two stages of approximation are required. The average experimental 
value for the particular conditions was 0-174, which indicates that only about 
60% of the photoelectrons enter the multiplying stages. 

(6) Effect of incident light pulse intensity. By using the neon lamp source and 
neutral filters the dependence of the half width of the pulse distribution on light 
pulse intensity was determined. Figure 2 shows in a striking pictorial way the 
appearance of statistical fluctuations as the light pulse intensity is reduced. The 
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Fig. 2. Multiplier output pulses for different light pulse intensities from neon lamp source. 
(V is pulse height, ¢ is time.) 

(a) 20000 photoelectrons per pulse; 7 

(6) 2000 photoelectrons per pulse; 7 

(c) 200 photoelectrons per pulse; n 

(d) 6 photoelectrons per pulse; n 


photographs give the appearance of the multiplier output pulses on a monitor 
cathode-ray tube for four different intensities. Figure 2(a) shows the uniform 
pulse height obtained when the number of photoelectrons produced is very large 
(~2 x 104) and for which number 7 is almost zero. For fig. 2 (d) the pulse intensity 
was so small that only about six photoelectrons were produced per pulse, giving 
a half width of 7=0-94. ‘The random thermionic emission pulses, due to single 
electrons, can be seen among the regularly spaced light pulses in this photograph. 
Figure 3 shows the linear relation found between 7 and the inverse mean light 
pulse intensity, given by the transmission coefficient of the filters which is 
expected from eqn. (3). 


(iv) Study of Pulse Spread due to the Luminescent Crystal 


Various crystals were used in these investigations, some being optically 
perfect and others having marked flaws and irregularities. ‘The half width and 
mean pulse magnitude were measured for standard photomultiplier conditions, 
using a 5 Mev alpha-particle source as excitor. After each crystal had been 
studied the neon lamp light pulse source was swung into position over the 
multiplier cathode and the values of half width and mean pulse magnitude 
determined. By means of this check any variations in the multiplier behaviour 
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become apparent. The values obtained for the neon lamp and for various — 
crystals are plotted graphically in fig. 4 to show the relation between 7? and 1/ Vo | 
It is seen that all values for the neon lamp and the perfect crystals lie on the _ 
same straight line which passes through the origin, that is, pulse spread is mainly _ 
determined by the photomultiplier characteristics. However, the points for 
imperfect crystals deviate considerably from this straight line depending on the 
degree of imperfection in the crystal. It should be noted that, as would be 
expected, all imperfect crystals give points lying above the straight line. 


° 10 20 3Q 49 so 60 70 


Fig. 3. Variation of half width 7 of output pulse distribution with light pulse intensity (given by _ 
transmission coefficient T of filters) of neon lamp source. 


Fig. 4. Variation of half width 7 of output pulses with mean pulse size V for scintillations from 
various crystals and for neon lamp light pulses. 


1. Neon lamp pulses; flawless crystals. 2. Scheelite. 3. 4. Anthracene in naphthalene. 
5. Stilbene in naphthalene. 6. Stilbene in naphthalene; crystals with flaws. 7. Thallium 
activated potassium iodide. 8. Copper activated zine sulphide. 9. Anthracene. 
10. Anthracene in naphthalene. 11. Stilbene. 


$4. CONCLUSION 


The above results are of importance when using the scintillation counter for 
energy measurements. ‘I‘hey indicate which factors are mainly responsible for 
the spread in pulse magnitude at particle energies in the range 1-10 Mev. These 
are the number of photoelectrons produced per scintillation and the optical 
quality of the phosphor crystal together with its intrinsic efficiency for converting 
particle energy to luminescence. ‘The effects of multiplier dynodes are seen to be 
relatively small compared with the above factors. 

It would seem very desirable to improve the photoelectron collection efficiency 
for the two types of photomultiplier investigated, and if possible to improve the 
intrinsic efficiency of phosphors. Since only 20% of the absorbed energy 
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appears as luminescence a fundamental investigation of the non-radiative 
processes by which the other 80% is lost would be useful. Some attempts to 
investigate the latter have been made recently by Schén (1951), by Broser and 
Warminsky (1951) for inorganic sulphide crystals and by Birks (1950) for 
anthracene. 
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The Current in the Electron Immersion Objective 


By L. JACOB 


George Holt Physics Laboratory, University of Liverpool 
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_ ABSTRACT. It is shown that for an electron immersion objective under the condition of 
constant cut-off voltage the emission current is uniquely defined by the cross-over potential. 
There is thus some theoretical foundation for the empirically observed law that the current 
at zero modulation voltage depends only on the three-halves power of the cut-off voltage. 


Sie UNC R OD UE OIN 

T is well known that the electron immersion objective, in common with all 
] other electronic devices which work under space charge limited conditions, 

obeys the three-halves power law. ‘The theoretical problem, however, of 
determining the importance of the various parameters, in connection with any 
desired modulation characteristic, is somewhat complicated. A contributing 
feature is the special geometrical construction, which is necessary in order to 
provide a rotationally symmetrical beam, of the current density required for 
practical operations. The main difficulty arises from the large values of the ratio 
of aperture diameter to cathode-aperture separation. ‘This may be as high as 
20 to 1, and thus emphasizes the ‘hole’ effect. For this reason it has been 
customary to base design data mainly on empirical results; these usually cover a 


| large number of variations carried out on model systems, either in demountable 


tubes or in separate processed units. In this way it was found that for systems 
with different geometry but with a common cut-off voltage V,, the current 
furnished by the system was approximately constant and depended only on the 
value of the cut-off voltage. 

Examination of the modulation characteristic revealed an empirical law of the 
form J = KV7?//2,, for the current at any voltage V, up from the cut-off point V,., 
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where the constant K has a value of about 3 ; this law appears valid to several per 
cent over a considerable portion of the curve. When V=V,, then [=KV,,%”, | 
which is the well-known three-halves power relation, for the cut-off voltage is 
directly proportional to the anode voltage. The purpose of the present paper is to 
seek an explanation for the constant value of current determined by a constant 
value of cut-off voltage. Different geometries were tried in which the cathode to 
modulator spacing was varied over a range of 10 to 1, keeping anode spacing, 
modulator aperture and thickness constant. ‘The system is shown in the figure. 


Cathode 1:8 mm diameter. 
| Modulator to anode 1 mm. 
Eo ore Cathode/modulator spacing—variable. 
Modulator aperture 1:0 mm. 


ie diameter 10 mm. 
e thickness 0:3 mm. 


Modulator 


Anode 


§2. SOME CROSS-OVER PROPERTIES 


It was shown earlier (Jacob 1950 a) that the axial potential distribution of 
the immersion system can be expressed as a sum of exponential terms involving 
certain constants of the system. Over the ranges explored this can be approxi- 
mately represented by the simple law V = Ae”, where A is a voltage scale factor 
connected with the form of the field, that is, it depends separately on anode 
voltage and modulator bias, while k depends on the geometry. On applying this 
to an analysis of the trajectory equation it was possible to derive certain properties 
of the cross-over (Jacob 1950 b). These may be summarized as follows : 


Distance of cathode’to cross-over =S8n/3Ry/3. ~~ = eee (1) 
Potential at cross-over = A expi(S7/34/3) 5 95 | ieee (2) 
Radius of cross-over =(2tan@/k) exp(—27/3+/3) ...... (3) 


where @ is the limiting angle of emission, within the paraxial range, of electrons 
leaving the centre of the cathode. Thisis taken to have aconstant value throughout 
the analysis. 

It may thus be seen that the potential at the cross-over simply depends on the 
value of the constant A. ‘This can be derived from the plots in the electrolytic 
tank for any set of field conditions. A typical set of values is given in table 1 fora 
range of variation in cathode to modulator, C/ M, spacing of 10/1. 


Table 1 


Modulator thickness 0-2 mm; modulator aperture diameter 1-0 mm; 
separation from anode 1 mm. 
C/M (mm). 0:05 0-10 0-30 0-50 
A (for 30 v cut-off at zero bias). 0-96 0-98 115 0-72 
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These observations indicate a tendency for the potential at the cross-over to 
remain constant over a substantial variation in spacing, provided the condition of 
constant cut-off voltage is maintained for each geometry. 

The variation for the larger spacings can be accounted for by an error of some 
10% in determining the cut-off voltage under the conditions in the tank, since the 
voltage scale had to be stepped down by a factor of 12. 

Independent confirmation of the constancy of the cross-over potential was 
obtained by actual ray tracing on a similar type of objective whose field plots were 
also determined in the tank (Einstein and Jacob 1948) for each variation in 
geometry. 

The actual position and size of the cross-over will, of course, vary from geo- 
metry to geometry, depending as it does on the k of the system according to (1) and 
(3) ; however, the gradient at the cross-over (whichis numerically k times its value of 
potential) does not vary greatly except for the closest spacing for which k has its 
highest value. As the potential at the cross-over is directly proportional to the 
potential on the anode, which in turn governs the cut-off voltage, we see that the 
cross-over voltage and cut-off voltage are proportional. The approximate 
relationship between anode and cut-off voltage can sometimes be calculated in 


# simple cases, e.g. where the anode and cathode are the end pieces of a flat pill box 


at potentials V, and V, respectively, spaced infinitely close to, but insulated from, 
the body which forms the modulating electrode kept at zero potential. The 
potential V at any point (2, 7) is given (Jacob 1951) in terms of the radius aand length 
2/ of the cylinder by 


J (kyr) cosh kyz J (Ror) cosh k,z ie 

k,aJ,(k,a) cosh k,/ ti k,aJ ka) cosh Rol 
J((Ryr) sinh kyz J (Ror) sinh ky 

+(Va-Vi)| po ea aer el + keiteayonk ent 


where J, and J, represent Bessel functions of zero and first order respectively, 
and k,, k, etc. are constants defined by the roots ofthefunction. For points onthe 
axis J, (k,,r)=1, hence the axial gradient at the cathode can be derived by differ- 
entiation of (4) with respect to z. The cut-off condition corresponds to zero 
gradient at the centre of the cathode, apart from initial velocity effects, i.e. the 
neutralization of the field which penetrates through from the anode, by the field 
from the modulator. In this case the zero equipotential just meets a central 
region of the emitting cathode, while the gradient over the remainder of the 
surface is negative. 
The relation between the cut-off and anode voltages for any system of such 
simple geometry is then given by 
Vi+V,  — X(cothk,,)/J,(k,.4) (5) 
VV, ~~ Sanh D(a) 


If V, is desired to have any value (in this case positive, and numerically equal to the 
cut-off value) the corresponding anode voltage V, can be found in terms of the 
geometrical factors (a, /), but these would require to be worked out for each 


geometry. 


V=(Vi+V2) 


§3. THE CURRENT IN THE BEAM 
In the following analysis it is assumed that the static distribution is not 
affected by the space charge, and the complicating factor of initial velocities is 
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neglected. We make use of the following well-known approximation for rotation- 
ally symmetrical fields: | 

r r* 
3 64 
and apply four fundamental conditions for space charge flow : (a) Poisson’s 
equation, (5) current density =space charge density x velocity, (c) v=(2eV/m)¥# 
and (d) use of relation V=Ae™. Take any element of space within the region of 
current flow, then 


VOLE Ay OE eee (6) 


loa OV o7V 
pa alin pa See oe ae a 
r or (- 7) Oz? gre (7) 
giving 
VEEP VPS = Fp Oe (8) 
and 
P= 77 REG er Et eee (9) 


This gives the volume density of charge at a point distant 7 from the axis in terms of | 
the potential on the axis and the k of the system. Now from (d) 
i { s (ie pany) eS {i (=)"a aie nay pe a (10) 

This is the current passing through unit area at the point (z,7). Now imagine a 
ring element of width dr to embrace this point : 

Current in ring element 7, = (277° k*/167)(2e/m)"(1 — 47? k?)12 V,3? dr, 

Current across section of beam j, = }(2e/m)#?V 3 i 2 73 k4 (1 —47? Raye dr, 
where the limit of 7 corresponds to the periphery of the beam. 

ik OF ice le 73k4(1 —7* k?/8) dr since r<k 
=C, V3? [47tk* —er' kG. 


Now, it is clear that V, can be taken anywhere along the beam axis, e.g. at the 
cross-over. In this case the limits of R correspond with the radius at this point, 
i.e. C, tan 0/k, where C, is a numerical constant given in (3). Then current over 
section at cross-over 


C,4tan*@ C,®tan®@ 
aot eciT Fog 8/2: | |e cg pe ae es 
= C1 Ver [ 4 48 | 


As C,, Cy, @ are constants on the simple theory, it is seen that the current in the 
beam depends primarily on the potential at the cross-over ; hence, for different 
systems which have the same value of cut-off potential, and hence of cross-over 
potential, the current should have a unique value over this range of geometries. 
This is found in practice. 

This relation should also hold when the beam is modulated. In this case it is 
necessary to determine the value of A corresponding to each modulation point. A 
typical set of values is given in table 2 (obtained from the tank plots) where M is 
the modulator voltage. 


Table 2 


M 0 —0°54 —1:0 —1°5 
A 0°43 0-30 0-21 0-04 
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There is, however, some uncertainty in the value of 0 when the modulation is applied 
because its effect is to decrease both the emitting area and the beam angle, and it is 
not clear as to how @ varies with modulation. This means that it is not yet 
possible to predict with any accuracy the theoretical shape of the full modulation 
curve. 
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ABSTRACT. At wavelengths outside the photographic range the determination of 
the optical constants of absorbing media involves the measurement of intensity changes 
on reflection at the surface of the medium. A review of previous methods for their 
determination is made and the experimental difficulties are indicated. A method is 
described which surmounts some of these difficulties by avoiding the direct measurement 
of reflection coefficients. The accuracy to be expected from the method and an example 
of its use are briefly discussed. 


§1. REFLECTION METHODS—GENERAL THEORY 

HE determination of refractive index m and absorption index k from 

measurements on reflected intensities has been described in general terms 

by Tousey (1939). In principle the method consists in measuring some 
characteristic of the reflection at various angles of incidence. ‘Then, if the 
measured values of the characteristic at angles of incidence 6,, @,, etc., for plane 
polarized light of azimuth ¢, be R,, R,... . etc. we have R= (ak, 1, $), 
R,=F(n, k, 0, 4), etc., where the function F is obtained from the appropriate 
Fresnel’s equation for the reflection. ‘Then the pair of values of m and k which 
can be obtained from the simultaneous solution of these equations are the 
q optical constants of the medium. Since the algebraic relations cannot be solved 
), directly for 7 and k the equations are solved graphically. 


§2. THE MEASUREMENT OF REFLECTION COEFFICIENTS 

Of the various characteristics of a given reflection, the one most generally 
chosen for measurement is the simple ‘reflection coefficient’, i.e. the relative 
intensity after reflection under given conditions of incidence and polarization. 
This method was described by Collins and Bock (1943) and has recently been 
discussed in detail by Simon (1951), who publishes some graphs relating the 
reflection coefficient to and k for angles of incidence 20° and 70°. The 
measurement of reflection coefficients, particularly at non-normal incidence, 
presents a number of difficulties which will be briefly discussed below. 


* Now at Telecommunications Research Establishment, Great Malvern, Worcs. 
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Perhaps the greatest of these problems are those associated with the large 
decrease in effective aperture of the reflecting surface as the angle of incidence — 
is increased. As in many cases where the optical constants of single crystals, etc., 
are being determined, the specimens available tend to be rather small, and since © 
in no case must the specimen itself become the limiting aperture of the optical 
system, it follows that the beam cross section used must be very small, with | 
consequent loss in intensity if beams of only small divergence are to be employed. 
This problem becomes particularly acute when the values of m and k are such 
that, to maintain the accuracy of the experiment, one reflection coefficient 
measurement must be made at 70°, or even higher, incidence. 

A second difficulty is that, for accurate reflection coefficient measurements, 
not only must the beams used have a constant intensity over their cross sections, 
but also the same portion of the photosensitive surface of the detector must be 
used at all angles, to avoid errors due to variation in photosensitivity across the 
detector surface. For accurate work (Kuhn and Wilson 1950) this involves the | 
use of diffusing screens in front of both source and detector, with consequent | 
further loss in intensity. 

Finally, it is necessary to arrange for the source intensity and detector 
equipment to be constant over the times (at least several minutes) required to. 
make measurements on a specimen at a given wavelength at, say, three angles 
of incidence and, in addition, to obtain a reading for the incident intensity or for 
the reflection from a reference mirror. 


§3. THE METHOD ADOPTED 

Most of these difficulties no longer exist if, instead of taking as the reflection 
characteristic the reflection coefficient for a given plane of polarization, one 
measures the ratio of the reflection coefficients for incident light polarized in 
and perpendicular to the plane of incidence. ‘The basic requirements of the 
optical system in this case are then simply (a) the polarizer must be such that the 
beam passing through it suffers no lateral shift on rotating the plane of 
polarization, and (6) if the light incident on the polarizer has not equal components 
in the two planes of polarization the ratio of the components must be known. 

If these conditions are satisfied, then even if the specimen is so small, or 
irregular in shape, that it becomes the limiting aperture at high angles of incidence, 
this introduces no errors since the aperture is the same for both planes of 
polarization. Motion of the reflected image over the detector surface as the 
incidence angle varies is of no importance since no absolute measurements of 
reflection coefficient need be made. Furthermore, the period of time for which 
the source and detector need be stable is only the time necessary to turn the 
polarizer through 90° and to take two meter readings. In addition, no absolute 
measurements of incident intensity are required. 

The formulae relating the ratio of the reflection coefficients of a substance for 
two planes of polarization to the optical constants can be derived from Fresnel’s 
equations. If the complex refractive index of the substance be N=n(1—ik), 
defined above, and if the ratio of the reflection coefficients be p?, then at angle 
of incidence 6 

o_ 4° +6? —2a sin 0 tan 0+ sin? 0 tan? @ 
Pp @ +P + 2a sin 0 tan 0+ sin? 0 tan? 0’ 
where N?—sin? 6 =(a—1b)?. 
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Curves relating p? to m and k for angles of incidence of 30°, 45°, 60°, 73° 
and 80° have been plotted from computations for values of n between 2-0 and 8-0 
and values of k between 0 and 2-0. The curves for 60° and 80° incidence are 
shown in figs. 1 and 2. From two or more measured values of p” at different 
angles of incidence, tables of values of m and k can be drawn up from these graphs 
and plotted against 7 and k as axes. In order to obtain well-defined intersections 
it is desirable to arrange the angles of incidence so that measurements of p” are 
made on both sides of the quasi-polarizing angle for the material. One such 
intersection for three measured values of p? is shown in figas: 


“0 01 02 03 0-4 05 06 


je 
Fig. 1. Curves relating p? to n and k for 60° incidence. 
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Fig. 2. Curves relating p* to 7 and Fig. 3. Curves showing intersection for three measured 
k for 80° incidence. values of p? plotted against n and & as axes. 


§4. TECHNIQUE, AND THE ERRORS ARISING IN THE METHOD 

The specimen is fastened to an adjustable mount at the centre of a hollow 
block whose temperature can be varied with coolants, etc. ‘This block is mounted 
on a shaft with a slow-motion drive and a vernier scale to measure angles of 
incidence to 1’ of arc. The detector is mounted on an arm rotating about the 
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same axis as the main shaft, and geared to it with a 2 : 1 ratio, so that, once adjusted, 
the reflected beam remains on the detector as the angle of incidence is varied. _ 
The whole is provided with a vacuum chamber for work with cooled specimens. — 

Measurements are at present being made at wavelengths in ‘the infra-red, _ 
and the polarizer employed is a pile of selenium films, constructed and used as 
described by Elliott, Ambrose and Temple (1948). 

The principal source of error in these measurements lies in the possible 
non-fulfilment of condition (a) referred to above, that is, in possible movement of 
the reflected beam over the detector surface as the polarizer is rotated. In this 
connection the selenium polarizer used is excellent, since the selenium films 
used are extremely thin, and the eight films used in this polarizer would only 

produce a lateral shift of the order of 254. For accurate measurements it 1s 
necessary for the efficiency of the polarizer to be high, especially when measure- 
ments are being made near the quasi-polarizing angle for substances whose 
absorption indices are low, since then the value of the reflection coefficient 
for the parallel component may be very small, and if any of the unwanted 
component is present in the incident light it will be strongly reflected and 
cause serious errors. 

As in the case of direct measurements of reflection coefficient, it is necessary. 
to make corrections to the observed values of p? if the incident light has been 
partially polarized by transmission through a prism monochromator. ‘This 
correction factor is determined by a preliminary experiment, using near-normal 
incidence reflection from a freshly evaporated aluminium surface, conditions 
under which the real value of p? is unity. 

With the system as described, some preliminary measurements have been 
made of the optical constants of lead sulphide and lead telluride. The values. 
obtained have been previously reported and discussed (Avery 1951). 

The values of m and k obtained are accurate to some 4°, arising from errors 
in the measurement of p? of some 1—2%. ‘These are the errors arising from 
the optical system: no estimate has been made of possible errors in 
n and k caused by methods of preparing the reflecting surface. A tolerance of 
+0-5° in the value of @ does not introduce comparable errors in the values of 
n and k, and this represents a valuable tolerance in the degree of flatness required 
of the specimen surfaces. 
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Primary Aberrations and Conjugate Change 
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ABSTRACT. Expressions are derived, in terms of the functions normally used in optical 
designing, for the variation with conjugates of the primary aberrations of any symmetrical 
system with isotropic object and image spaces. The possibilities of invariance of the 
aberrations, or their simultaneous correction for more than one magnification, are discussed. 


S12 IN CRODWCRLON 

HE proposition that “it is impossible, by means of any combinations of 

reflexions and refractions, to produce a perfect image of an object at two 

different distances, unless the instrument be a telescope” was given by 
Clerk Maxwell in 1858, and since that time various aspects of the laws governing 
the variations of aberrations with changes of magnification have been 
investigated by a number of writers, among them Herzberger (1931, 1946), 
Boegehold and Herzberger (1935), T. Smith (1948) and Grenat (1949). In most 
cases the investigation has been based on the expansion of the Hamiltonian 
characteristic function or Bruns’ eikonal, the mathematics not always being 
simple, and where analytic expressions for the variation of the primary 
aberrations have been deduced these have generally been expressed in terms 
of special functions not normally calculated in optical work. In this paper it is 
shown that the variations of the primary aberrations with conjugates may be 
deduced from extremely simple considerations for any systems having axial 
symmetry and isotropic object and image spaces, the result being given in terms 
of the normal aberration functions. ‘These expressions readily give the various 
possibilities of, and conditions for, invariance of one or more of the aberration 
coefficients, or their simultaneous correction for two or more magnifications, 
and the physical realization of these theoretical possibilities is also considered. 


Se OR 

The first order aberration coefficients for a symmetrical optical system of 
spherical reflecting or refracting surfaces are simply related, the successive 
coefficients for the oblique aberrations Sj; (coma), Sj; (astigmatism), 
Sy (distortion), being derived from Sj, the coefficient relating to spherical 
aberration, by successive multiplication by a factor A/A (defined below), with 
the addition of the Petzval field curvature in the last case; for an aspherical 
surface, the additional oblique aberrations are similarly formed by successive 
multiplication of the axial aberration term by a factor HE. The total first-order 
aberration of the system is the sum of these surface coefficients. Longhurst 
(1950) has shown that the coefficients relating to spherical aberration of principal 
rays for a system of spherical surfaces may be expressed in a form representing 
an addition to the above series, and by an obvious extension the corresponding 
coefficient for aspheric surfaces, denoted below by Sy;, may be similarly 
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included. This coefficient is contained in the expansion of the eikonal given by 
Schwarzschild (1905), but is disregarded by him in his discussion of the primary 
aberrations. Using a notation similar to that given by Hopkins (1950), the 
first-order aberration coefficients for a single surface may be written 


S, = A*yA(u/n) + 8aytA(n) =0,+0, 


3 
Sy Loe +o,HE, Sy =9, (5) + PHOS +o,(HE)? 
Se (4) +0,(HE), Sy =HE(Sy+PHA5) 

3 
where y is the intersection height on the surface for a paraxial ray at the conjugates 
considered, wu its angle with the lens axis, 7 its angle of incidence, m the refractive 
index, P the Petzval field curvature, H the Smith—Helmholtz invariant, 
A denotes the change on refraction, A =ni, a is the coefficient of fourth power 
asphericity, given by x=ay*, where x is the departure of the surface from 
sphericity ; a bar over any symbol indicates that it relates to a paraxial principal 
ray through the stop centre instead of to the paraxial ray from the object. 

Again following Hopkins (1950), A/A =ni/ni=H(E+1/Ay) and HE=jy/y, 
the quantity E so defined being generally derived in practice from data of a 
paraxial ray from the-object, using Seidel’s difference formula. 

Clearly the aberrations of pupil imagery, for which the object plane (field 
stop) is the effective stop, can be expressed by equations identical with (1), with 
those quantities relating to ray positions with and without bars interchanged, 
n, a, and P remaining the same in each case; and by definition HE=1/HE, 
ip Pee ul Denoting these first-order aberration coefhicients of pupil imagery 
by S,, Sy etc., then by definition S,=Sy, Sy,=S, and by simple algebra it 
may be shown that 


Sy—Sy=HA(u), Sp—Syr=HA(ua), Sij—Sy=HA(@). ...... (2) 
For example, from eqns. (1) 
Syi= 5; = HE(Sy + PHA}), 
and since Sy =S,HE+ HAA(a/n) 
and S,= A’yA(a/n) + 8aytA(n), 
Sy= Sit AH { Py — A(a/n)} 


Sy—S,=HA | ni(- a *) | = —HA(#@—7) = — HAR, 


n nN 


and similarly in the other cases. | 
On summing through the complete system the intermediate terms on the | 
right-hand side of each of eqns. ve cancel, giving 


a Sct Sy =A(u?—u,?) 
S Sire RTE tn POI (3) 
YS —LSy =H(a?—4a,2), 


where wu, uj, %, m3, are the angles with the axis of. the paraxial ray and paraxial | 
principal ray in the object and image spaces, ‘These equations have applications 
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in other connections than the present one; thus the third equation shows how 
the coma of stop imagery S};, which is of importance in considerations of the 
intensity of illumination in oblique imagery (Slussareff 1941), depends only on. 
the primary image distortion, the angle of field and the pupil magnification. 

Now the change in the primary aberration coefficients produced by change 
of stop position is very simply expressed, since the effect is to add to HE, in 
eqns. (1), a term constant throughout the system. A change of stop position for 
pupil imagery is equivalent to a change of object position, so that, by applying 
eqns. (2) to the equations showing the variation of S,,, Sy; and Sy with object 
position, the changes of Sj, Sy; and Sy with conjugates are readily derived. 
If the object be considered to be moved such a distance as to add a term HdE 
to HE in the expression for S},, Sj, and Sy, then, from the equations for these 
coefficients analogous to eqns. (1), the new values S,,*, Sy;*, etc., are given, 
after some reduction, by 


ESyy* =US pq, + 2ZHSE DS, + (HSE)? ZS, 

USy* =USy+ HSEX(3S,,,+ PH?) + (HSE)223S;,+(HSEP=S, 

BiSy* =D Sy, + HOE X(3S_+ Sy) + (HSE)? (3 Si +35 + 2PH?) 
+ (HSE)? X(Sy+3Sq) +(HSE)tES;. 


BS,* =35;,+ BESS, | 
: 


Equations equivalent to these have been given by a number of writers; they 
are given in this form for an aspheric system by Wynne (1949). 


Eigostcas 


If A (fig. 1) be the axial object point for the original imagery, and P the pupil 
point in the object space, the paraxial ray from A, AB, making an angle w with 
the axis, and the paraxial principal ray DP making an angle # with the axis, then 
if the object plane be moved to F, the plane through F normal to the axis 
intersecting AB in G and DP in E, the angle BFP will then be u*, and the 
value of HSE corresponding to this change of conjugates (keeping H constant, 
i.e. for a constant pupil size and field angle) is given, by definition, to the 
order of approximation required, by 


HSE=GF/FE=(u*¥—u)/i. 2 neeee (5) 
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Substituting Sj,, Sy and Sy from eqns. (2) for Se, Siz and Sip» and u* from 
eqn. (5) in eqns. (4), the change of the first-order aberrations for a change of | 
conjugates are given by 


USy*' —USy =HSEXUS, 
DSin* — VS yy, = HSE X(Sy + Sy) + (HSE)*=S, 
= HSEX(2Sy + HA®) +(HSE)2=S, 
Sq* —ESp =H D2S n+ Sun PH) + (ASEY E(Sy +25q) 
+ (HSE) ZS, 
= H8EX(3S yy, + PH? + HAui) 
+ (HSE)? X(3 Sy + 2HA@) + (HSE)3 US, 
US,* —US, =HSEX(3S y+ Sy) + (ASE)? X(3S_q7 + 3Spqy + 2PH?) 
+ (HSE)? 3(Sy+3Sq) + (HSE) ZS, 
= HSEX(4S,,+ HAu?) 
+ (HSE)? X(6Sy;;+ 2PH? + 3H Aun) 
+ (H8E)3X(4Sy+3HA@)+(HSEYIS;, . | 
The proof of eqns. (3) and: (6) given above is valid for any symmetrical 


optical system of refracting or reflecting surfaces. The same equations may 
be shown to be generally valid for any symmetrical imaging system having the 


) 


object and image space isotropic. Considering first the case of axial spherical 
aberration, let A, A’ be the original paraxial object and image points and P, P’ the 
paraxial entrance and exit pupil points, all lying on the axis APP’A’ (fig. 2), 
and AB, B’A’ be the incident and emergent paraxial rays. If the axial aberration 
is required of a point F, whose paraxial image is F’, then drawing BF to meet 
the plane normal to the axis through A in C and taking circles through F, centred 
on P and C, to intersect CP in Q and R and similarly in the image space; denoting 
optical path lengths by square brackets, and the aberration associated with the 
imagery of the point A by wy, and similarly for the other points, then the 
aberration required, wy is given by 


wyp= —[FB...B'F ]+[FP... PF] and w,=—[FB... BF] >( RP Pe 

ip = - [OP ok OF eR ee Pee | 

or Wy =@y+wp— An(QR). 
If now the first-order values be substituted in eqn. (7), the aberrations for 

the plane AC being assumed known for a field angle #, then if the angle APC =y, 
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and again following Hopkins’ formulation for the first-order aberrations of a 
symmetrical system of spherical surfaces, which is also the form of first-order 
aberration for any symmetrical system, 


wp = 32S {(y/t)! 

wo = 28S p+ $Spry/h+ 33S q+ PH?)(y/d)? + $Sy(y/a)3]. 
From fig. 2, if the angle BFP =u*, QR=AF sec u* + FP— AP sec y. From the 
geometry of the figure, FP tanu*=AP tanu, AF tanu*=APtany. If H be 
the Smith—Helmholtz invariant for the initial aperture angle u and field angle a, 
then 7A Puy = Hy/a, the corresponding function for the field angle y, is also 


invariant. Substituting in the equation for QR and expanding the trigonometrical 
functions, to the order of approximation required, 


—AnQR=3HE | w+3 * uii+3 (1) | ‘ 
Substitution in eqn. (7) then gives 
ES,* ==S;4+ (y/#) D[4.S + Aw] + (y/a)? N[6.S,, + 2PH? + 3HAua] 
+ (y/i)? S[4Sy + 3HA®] + (y/a)4 US). | 


This latter equation is identical with the last of eqns. (6) since HSE =y/da. 
Alternatively, the aberration of the point P could be expressed in terms of 
@,p and w,, instead of we and w», an equation analogous to (7) resulting. By 


Fig. 3. 


inserting the first-order aberrations in this case, and_ equating coefficients in 
this and eqn. (8), the relationships between S;; and Sy etc. given in eqns. (2) 
are obtained, so that these too are valid for any symmetrical system having 
isotropic object and image spaces. 

The first-order oblique aberrations for a second pair of conjugates in terms 
of those of the first is given by a similar procedure. ‘The simpler case of 
aberrations in a meridian plane will be considered here, though the same method 
is obviously applicable to the general case. Let C, C’ in fig. 3 be paraxial object 
and image points on the initial conjugate planes PP’, BB’, the paraxial pupil 
centres and edges, respectively, and let F, lying on the line CP, be a point on a 
new object plane; then joining FB to intersect the original object plane AC in D, 
and adopting a similar procedure to that given for the case of axial aberration, 
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the aberrations of the imagery of F may be given in terms of wp, the 
aberrations wpq and wpg,,) of the axial pupil-point P for aperture angles — 
APC=a and APD=(a#+y), and the quantity Am(QR) where Q and R are the 
intersections with DP of circles through F centred on P and D, 


Wy = Wp + Wpg 4) — @p> — An(QR). ccc (9) 


Inserting the first-order values and equating coefficients, eqns. (6) result. 

The variations of the first-order chromatic aberrations may be treated 
similarly, and since in this case the order of magnitude considered is that of the 
second power only of the aperture and field, the quantities An(QR) are 
negligibly small. Denoting the chromatic aberration of the imagery of C for 
some specified wavelength variation by w,=4$UC,+UCyq where 42C; is the 
longitudinal and UC), the transverse chromatic aberration, and writing $2Cy 
for the longitudinal aberration of the pupil point for a field angle a, then by 
eqn. (9), to the order of accuracy required, 


= tg 5 5 
2UCy* + BCy* SEC eg = +22C) [ee or 1| +3HA Ge 


u 


or, equating the parts involving equal powers of the field angle, 


: 5 ote 
UCy* = UCy + (y/W)UC, 


where, for a symmetrical system of spherical surfaces (see, for example, 
Hopkins 1950) 2C;=XAyA(dn/n), UCy==UC,A/A, dn being the change of 
index of each medium for the wavelength range considered, and as before, in 
eqn. (10) y/a#= HSE. 

The variation of the first-order aberrations with stop position may readily 
be treated by the same general method, and the first three of eqns. (4) may be 
shown in this way to be valid, not only for any symmetrical system of reflecting 
or refracting surfaces, but for any system having axial symmetry and isotropic 
object and image spaces. 


§3. CONCLUSIONS 
(i) Invariance of Aberrations 


The conditions for the invariance of the Seidel aberrations for changes of 
conjugates follow immediately from eqns. (6) and, while each aberration 
separately may be made zero and invariant for all magnifications, it is clear 
that this cannot be achieved for all the coefficients simultaneously in a lens of 
finite power. Various combinations of invariance are possible. 

Distortion is invariant if the spherical aberration of the pupil point is zero, 
and astigmatism if in addition (Sy+HA#)=0; a lens having invariant 
astigmatism has therefore invariant distortion, and this astigmatism can only 
be zero if the stop is in the nodal plane, so that Ad? =0; this problem was partially 
treated by Merté (1920). For a lens with its stop at the node, HAau may be 
written —H?/nf’, where n is the refractive index of the object space and 
f’ the image space focal length. For a system with object and image space 
of air the three conditions for invariant coma are then S,=0; 3.S,+HA#=0; 
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3S),+ PH®?—KH?=0. For a system having zero and invariant distortion and 
astigmatism the first two conditions are satisfied, and the third, which is the 
condition for the coma to have a stationary value, reduces to equality of power 
and Petzval curvature, and cannot be achieved for a lens of finite power having 
a flat field. 

The four conditions for invariance of spherical aberration are S,=0; 
4S\+3HA@=0; 6Spy;+2PH?+3HAui=0; 4Sy+HAu2=0. Since the 
spherical aberration is independent of the stop position, the aberration of the 
pupil point is necessarily included in the range of invariance, and hence the 
coefficient S; can only be invariant if it is zero. Since Av? varies with the 
conjugates, being zero only at unit magnification, the last condition requires that 
the coma can only be zero at this magnification, corresponding to the 
incompatibility of the sine condition and Herschel’s condition at other 
positions. : 

The combinations of invariance of the primary aberrations which are 
theoretically possible are also physically realizable. Thus process lenses are 
made having invariant and corrected distortion astigmatism and field curvature ; 
such lenses have coma changing linearly with HSE, and spherical aberration 
stationary at unit magnification. For thin lens systems at the stop, distortion 
and astigmatism are invariant, and the number of conditions for invariance of 
coma and spherical aberration reduces to one and two respectively. ‘The 
condition for invariance of coma is satisfied, for example, by a thin Mangin mirror 
system of refractive index 1-620 whose air—glass surface is in the aplanatic 
condition at unit magnification, where the spherical aberration is corrected, 
and stationary. The conditions for invariance of spherical aberration may be 
satisfied by giving a thin system of mirror and meniscus lens of appropriate 
powers to make 6S,,,;+2PH?—3KH zero a shape to correct the spherical 
aberration of the whole. 

For the primary chromatic aberrations the invariance conditions given by 
eqns. (10) are very simple. ‘Thus transverse chromatic aberration is invariant 
if the chromatic aberration of the pupil point is corrected; and if the object and 
image spaces are of the same material, longitudinal colour is stationary if the 
transverse chromatic aberration be corrected, and invariant if, in addition, this 
transverse aberration is invariant. 


(ii) Correction for more than one Magnification 


Since spherical aberration is normally corrected for four sets of conjugates, 
coma for three, and astigmatism for two, the conditions that these conjugate 
surfaces may coincide for different aberrations may be formulated. Since 
distortion is normally only corrected for one pair of conjugate planes and, if it 
is invariant, the number of corrected magnifications for each of the other 
aberrations is reduced by one, distortion correction is not compatible with the 
conditions to be considered. 

Considering first the simultaneous correction for two magnifications of 
astigmatism, coma, and spherical aberration, since, when these are corrected, 
they and Sy are all invariant with stop position, the stop may, without loss of 
generality, be taken at the node, S; then referring to the aberration at this stop 
position, and the conditions required, from eqns. (6), are then 


2Sy2 = S;( PH? — H?/nf'),2SyH?/nf’ =—Sy;HAw. ...... (11) 
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It follows that the second conjugates for which the lens is corrected must 
have an angular magnification which is the reciprocal of the first, and that the 
distortion is equal and opposite in the two positions, and that the coefficient S; 
must be of opposite sign to the power of the system if the field curvature is zero. 
For a lens corrected at a given angular magnification, the values of Sy and S; 
required to give correction at the reciprocal magnification may thus be calculated 
and, for the practically interesting case of a lens in air corrected for field curvature 
and a real object and image, the resulting values are so large as to be unattainable 
in systems of normal construction, requiring surfaces which are either extremely 
deeply curved or at relatively large distances from the nodes; thus, Sy and S; 
have their minimum values for real conjugates if the lens is corrected for an 
object at infinity and at the front principal focus, and in this case Sy = + H?/y”, 
S; =2H?*/y*K, corresponding, for a field angle 7 =, to a distortion of 87:5°% and 
a spherical aberration of the pupil point of 600 wavelengths per inch of 
focal length. 

In considering this problem of a system giving aberration-free imagery for 
two pairs of conjugates, Hertzberger and 'T’. Smith, in papers already referred to, 
have stated that Maxwell’s Proposition [X was in error; this proposition is 
that “It is impossible, by means of any combination of reflexions and 
refractions, to produce a perfect image of an object at two distances, unless the 
instrument be a telescope, and /=n=yp,/y.,m=L”’ (l,m, and n being the transverse, 
longitudinal and angular magnifications, and 14, 2 the refractive indices of the 
object and image spaces). ‘This disagreement appears to be due to the fact 
that Maxwell’s definition of a perfect image, given at the beginning of his paper, 
is one having perfect definition and freedom from field curvature and distortion, 
which is not the case for the systems considered by Herzberger and Smith. 
With Maxwell’s definitions, his proposition is correct. This does not of course 
invalidate Smith’s refutation of Maxwell’s later statement that “‘we cannot get 
perfectly focused images of an object in two positions even at the expense of 
curvature and distortion”’. 

A further case of possible interest is the simultaneous correction of coma 
and spherical aberration for three sets of conjugates. Here again, if S; and Sy, 
are zero they and Sj,;; are unchanged by a change of stop position, and this may 
therefore be taken at the node without loss of generality, Sy and S; then 
referring to the values of these coefficients for this stop position, and the two 
conditions to be fulfilled are 


S,(3 Sty + PH? —2H?/nf’) =3Sy, Sy(3Sq_-+ PH? — H?/nf')=S,;HAuw. ...(12) 


If Sj;,;=0, these reduce to equations which are clearly incompatible with 
eqns. (11) for a system of finite power and flat field, so that such a system cannot 
simultaneously have a flat field anastigmatism correction for two pairs of 
conjugates and be aplanatic for a third pair. It may be shown that, for real — 
object and image distances for all three conjugates, Sy and S, have minimum 
values if for one set of conjugates the object be at infinity and the other two 
coincide at the front principal focus, or vice versa, and in this case Sy has the same 
value as for the anastigmat corrected for two pairs of conjugates, and S; has 
three-quarters of the value for that case; as before, therefore, these are so large 
as to be virtually unattainable practically. 
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Propagation of Ultrasonic Waves in Vapours near 
the Critical Point 
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ABSTRACT. An experimental study is made of the velocity and absorption of ultrasonic 
Waves in carbon dioxide and ethylene at frequencies between $ and’2 Mc/s and at pressures 
up to 100 atmospheres in a variable-path acoustic interferometer. At pressures above 
10 atmospheres the results are independent of frequency (in this range) in the gaseous phase. 
The velocity and absorption coefficient reach a minimum near the condensation point, but 
especially at the critical point itself have high and indeterminate values when the change of 
phase takes place. The ultrasonic viscosities calculated from these measurements are of the 
order of 1000 times those derived from transpiration or low-frequency oscillation experiments. 


Sel © Dei Ga TON 

XPERIMENTAL work on the propagation of ultrasonic waves in various 
gases has given results for ultrasonic velocities and absorption 
coefficients which cannot be adequately explained by the ‘classical’ 
theories of sound propagation. Generally, for given pressure and température 
conditions, there is, at certain frequencies, a dispersion of velocity accompanied 
by a sharp increase in absorption at that frequency. Experiment shows that an 
increase in pressure or temperature displaces the dispersive region to higher 
frequencies. However, with few exceptions, previous workers on the effect 
of pressure on ultrasonic velocity and absorption in vapours have covered only 

a limited range of pressure about the normal. 

It was the intention in this research to make measurements in one or two 
vapours at such temperatures and pressures that the transition to the liquid state 
would be covered and, in particular, to investigate the variation in the absorption 
coefficient « in this region. It is well known that the coefficients of variation of 
viscosity with temperature run conversely in the two states, so that information 
of this kind should show whether the viscosity which characterizes ultrasonic 


* Now at the University of Liverpool. 
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absorption follows a similar trend. For convenience, carbon dioxide and ethylene 
were chosen as they have critical temperatures not far removed from the 
temperature of a laboratory. The research was conducted as a study of 
ultrasonic propagation in the neighbourhood of the critical state and not, 
intentionally, as a. study of the molecular properties of individual vapours. 

Most of the experiments were made on carbon dioxide, with a small number _ 
on ethylene. ‘The well-known dispersion of velocity V in carbon dioxide at 
atmospheric pressure was absent at higher pressures over the range of frequencies 
used here. At all three frequencies, single isothermals for velocity plotted 
against the quotient frequency/pressure were obtained. Ethylene also shows 
some dispersion at pressures near atmospheric, but to a much less extent. In- 
fact, both the vapours used in these experiments showed similar ultrasonic 
properties, having regard to the location of their respective critical points, so 
that extensive work on ethylene proved unnecessary for our purpose. 

Previous work on these vapours under high pressure is confined to velocity 
and a few frequencies only. Spakowski (1935) made measurements in carbon” 
dioxide near the critical point at a frequency not stated, although from the 
description and size of the apparatus one judges it to have been in the audible 
range. Herget (1940) has measured velocities at one ultrasonic frequency of 
about 280 kc/s in carbon dioxide (5 to 98 atm) and ethylene (35 to 75 atm) at 
temperatures near their respective critical points. Finally, Noury (1946) quotes 
some velocities in what he describes as impure ‘war-time’ carbon dioxide. 


Fig. 1. Section of interferometer. 


§2., APPARATUS : 

A single quartz acoustic interferometer was used for the measurement of 
velocity and absorption in gases and liquids under high pressures. 

Figure 1 shows the final design of pressure vessel used. The main vessel 
contains the quartz transducer A, facing the reflector on its piston B, in a 
chamber to which the gas is introduced by a port C. The piston is moved by 
the large screw-head D whose setting is measured optically. There is a second 
gas chamber and piston of suitable size to compensate for changes in volume 
as the main piston moves. The whole apparatus is surrounded by a lagged 
enclosure. 
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A platform into which the crystal holder is locked is mounted in the vessel 
by three spring-loaded bolts. This permits parallel adjustment of the transducer 
face and reflector. 

The sources were X-cut crystals with plated electrodes, mounted so that the 
quartz merely rests in a tufnol cup ona brass bed-plate machined as flat as possible, 
heavily chromium-plated and then polished to a near optical flat. A pigtail of 
0-001 in. diameter platinum wire makes contact with the upper surface. 

The crystals were discs cut wide relative to their thickness to limit errors due 
to diffraction. he ratio of radius R to thickness (or wavelength) did, in fact, 
increase at the highest frequencies. Losses due to diffraction, if they had occurred, 
would have been apparent as an increase of « at the lower frequencies, but such 
was not the case. 

Though an examination of the wave-forms during a ‘run’ will show whether 
or not source and reflector are accurately parallel, an initial test was made in which 
this system formed part of a Michelson interferometer so that the parallel 
movement would be assured to within a light wavelength. ‘This is unnecessarily 
precise for the present purpose, but the necessity for parallelism cannot be 
over-emphasized if accurate values of absorption are to be attained, since 
otherwise satellite peaks or cols in the reaction wave-form are obtained. The 
initial adjustment to ensure that the reflector face is parallel to the quartz (before 
the base-plate is clamped on) is made with a depth gauge and the final test by 
observing the acoustic wave-form. It may not be possible entirely to eliminate 
satellites by such parallel adjustment if transverse resonances of the fluid in the 
tube are excited, as Bell (1950) has shown, but the fact that the accepted values of 
absorption in air were obtained in this apparatus indicated the absence both of 
this effect and of variations in phase and attenuation coefficient due to the 
consequent three-dimensional standing waves. 

The crystals were driven by a Hartley oscillator of conventional pattern, 
which gave stability within +1 c/s over the time of one run (about 7 minutes) 
| up to a frequency of 4 Mc/s. A frequency meter, checked from time to time 

against B.B.C. transmissions, accurate to 1 in 10’, was loosely coupled to the 
© oscillator. The frequency changed slightly under pressure. Consequently 
the frequency ascribed to each isothermal curve of V and « in the sequel is to be 
considered as the mean for the curve in question. 


§3. FILLING THE VESSEL 

Gas in cylinders supplied with special attention to the absence of air as an 
adulterant is stated to be better than 99-7°% pure. As explained in the Introduc- 
tion, our concern was to have the gas free from any adulterant which might 
cause dispersion of velocity specific to its presence and consequent 
‘relaxational absorption’. 

The possible impurity most likely to do this is water vapour. However, 
attempts to remove the last traces of water vapour by drying agents or to add 
further water vapour (in one special test) produced no change at s.t.p. in the 
dispersion or absorption over the working range of frequency. At higher 
pressures, as will appear, these phase and attenuation factors become independent 
of frequency, at least as far as the gaseous state is concerned, in the range of 
frequency employed. Nevertheless, these s.t.p. tests—by comparison with the 
known results of other workers in carbon dioxide—serve as useful checks on the 


purity of the gas supplies. 
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The bomb was evacuated and connected to the cylinder. To obtain 
pressures up to 120 atm the gas cylinder was carefully heated before it was 
cut off from the bomb. The Bourdon gauge, duly calibrated, gave the pressure 
of the gas inside the bomb. 

‘Temperatures over the range 0°—40° c were attained in the lagged enclosure 
by circulating air which had been fanned over electric heaters (for carbon 
dioxide) or drikold (for ethylene) as the case may be. ‘The thermostat for the 
enclosure is a toluene regulator fitted with a Sunvic proportioning head. The 
regulator controlled the heater through an electronic relay. With this 
arrangement it was possible to hold the gas temperature, shown by thermo- 
couples, inside the bomb to 0:02°c. A high grade of stability in pressure and 
temperature is, of course, necessary near the critical point. 


§4. METHOD OF MEASUREMENT 

The conditions under which the attenuation in the gas column can be derived 
from the logarithmic decrement of the peaks of anode current J in the driving 
circuit as the reflector is moved have been laid down by Hardy (1943). The 
somewhat different conditions for a liquid—due to the more nearly matched 
impedance of fluid and solid reflectors—have been discussed by Fox (1937). _ 
Hardy shows that errors in the absorption coefficient « will not exceed 2% if 
the ratio 5//I is less than 0-03, where 85J is the change in anode current as the 


Fig. 2. Peaks in anode current with change of path length. 


reflector is moved from a position of maximum to a position of minimum. To 
measure small changes in J a shunted galvanometer with ‘ back-bias ’ arrangement 
was used. Battery supplies were used throughout. 

Automatic recording of the galvanometer on a moving-paper camera was used 
to eliminate personal errors in reading the moving galvanometer light spot and 
to permit a constant check of wave-form. 

The cylindrical lens in front of the camera has a series of lines engraved on its’ 
face. With a suitably disposed small filament lamp the photographic paper can 
be fogged except where these engraved lines cast a shadow. By arranging that 
this light fogs the paper before the exposure due to the galvanometer light—a 
guestion of orientation of lamps only—the initial photographic threshold of the 
paper is exceeded and an adequate exposure with the available reflected light 
obtained. Variable resistances and switches in series with these two lights 
give control of exposure times. 
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Figure 2 gives a typical run as recorded by this means. The speed of 
displacement of the reflector was considerably reduced at the beginning, middle 
and end of the trace. (The slightly smaller peaks obtained at other points in the 
trace are due to the motor ‘overrunning’ the response of the galvanometer.) 

The absorption coefficient between two peaks corresponding to two reflector 
distances /, and /, was calculated from 


df, _ sinh 2a(/,+8) _ 
dl, sinh 2a(/,+8) _ 


for sufficiently large values of «/ (8 depends on the reflection coefficient). 

In the actual mode of operation of the apparatus the gas was introduced 
at the highest pressure, left to attain the desired temperature, the various circuits 
put into operation and a run made at this pressure. ‘Then the pressure was 
reduced by allowing some gas to leak slowly into an evacuated vessel and another 
run made, when pressure and temperature stability was regained. 


exp 2a(/,—1,) 


§5. RESULTS 


The velocity results in carbon dioxide are plotted in two stages. In 
fig. 3 are shown the (V, p) relationships: straight lines well below the critical 
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Fig. 3. (V,p) isothermals in carbon dioxide. 


point and curving down to cusps at the critical point itself. ‘They are in general 
agreement with Herget’s results at 280 kc/s. ‘The frequencies used were $, | and 
2 Mc/s; the velocities were found independent of frequency in the gaseous 
phase at all pressures above 10 atmospheres. 
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The absorption results in carbon dioxide are shown in fig. 4 andinethylene | 
in fig. 5 in terms of the absorption per wavelength ». As already stated, fewer _ 
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measurements were made with ethylene, but fig. 5 shows typical results for both 


V and » in this substance at $ Mc/s. 


The velocity in this gas is also found — 


independent of frequency, when Herget’s few results are added to ours. 
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Fig. 6. («/f*, p) curves in carbon 


dioxide at 26:5° c. 


Figure 6 shows the constancy of «/f? in the gaseous régime, but discrete curves 
for the different frequencies in the liquid state, in which, however, the accuracy 
of measurement is less than in the gaseous state. 
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When by chance the conditions corresponded exactly to the point of 
transition, the reaction on the source failed completely, due to the very high 
absorption; this may be ascribed to scattering by the striated masses of 
indeterminate fluid known to exist in this regime. Proceeding, as usual, from 
the liquid towards the gaseous phase, sometimes a rise, sometimes a fall was 
encountered in just before the transition (cf. inset to fig. 4). In the former case, 
the absorption after reaching a value impossible to measure dropped sharply 
to a point on the curve proper to the gaseous phase. Schneider (1951) has drawn 
attention to the very large attenuation at the critical point itself. 


§6. COMPARISON OF ULTRASONIC AND STEADY-FLOW VISCOSITIES 

According to the Helmholtz—Kirchhoff formula, the absorption coefficient 
% in a gas at pulsatance w should be related to the kinematic viscosity v, thermal 
diffusivity « and ratio of specific heats y by 


When, however, we attempt to deduce from our data the variation of v with 
p using this formula we find very large values (fig. 7) in comparison with values 
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of v deduced from transpiration or low-frequency oscillation values. As the 
state of the substance was liquid or near-liquid, the second term in the brackets 
was neglected in comparison with the first in calculating these values of v. 
Phillips (1912) gives values of viscosity in carbon dioxide under pressures 
and temperatures close to the critical point. ‘These were obtained from flow 
experiments in capillary tubes. As there seemed to be no corresponding 
published values for ethylene (except at constant density) it was decided to obtain 
them from the damping of a vessel filled with the fluid oscillating about an axis 
of suspension in a vacuum. ‘This method, suggested by Helmholtz and 
Piotrowski (1860), was refined into an absolute method by Andrade and Chiong 
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(1936). A comparative method was used, with the gas cylinder itself (10 in. x 3in.) 
held by a bifilar suspension so that it could be gently set in oscillation—through the 
intervention of electromagnets—about its axis, set vertically. ‘The whole was 


mounted in an evacuated glass cylinder and this again in a thermostat, the 


damping of the oscillations being recorded optically. Afterwards, the cylinder 
was filled with a fluid of known viscosity and the damping again measured for the 
purpose of calibration. 

The results of such experiments on ethylene and carbon dioxide at the critical 
temperature are shown in fig. 8, using Amagat’s density values. It is at once 
apparent on comparing these curves with those of fig. 7 that the viscosities 
deduced from the ultrasonic attenuation using Kirchhoff’s formula (1868) in 
this region are about 1000 times greater, though the general trend of the curves is 
reproduced. 

It is interesting to note that although the actual values of the coefficient for 
gases under pressure are much higher than those deduced from Kirchhoff’s 
theory they follow the same variation with pressure, as is indicated by his equation. 

It may be remarked that ‘condensed gases’ at low temperatures such as liquid 
argon and helium show absorption coefficients of the same orderas do these vapours. 
in the critical region. 
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ABSTRACT. The contact of electrolytes has a marked effect on the mechanical behaviour 
of single crystals of cadmium and zinc which have been exposed to the air. This is due 
to the disruption, or formation, of a thin surface film, probably of oxide or hydroxide. 
If ordinary single crystal cadmium wires are loaded so as to flow very slowly, the contact 
of solutions containing free cadmium ions increases the rate of flow : a similar effect exists 
in the case of zinc treated with solutions of zinc salts. If the surface of the wire is cleaned 
by preliminary thermal evaporation there is no such effect, establishing that it is due to 
a surface film normally present. Solutions of cadmium or zinc nitrate stop the flow and 
uf raise the critical shear stress appreciably, which is attributed to the formation of a surface 
film, probably of hydroxide. With the cadmium nitrate there is a temporary increase 
_} of flow before the cessation, which is attributed to the cadmium ion, for there is no such 
© effect with the zinc nitrate. With polycrystalline wires there is no electrolyte effect. The 
work described supports the view that a thin oxide or hydroxide film increases the 
* mechanical resistance of single crystal wires, and that with cadmium and zine wires such 
| a film is normally present. 


$1. INTRODUCTION 

HAT the mechanical properties of metal single crystals are dependent upon 
the condition of their surface was established by Roscoe (1936) with single 
: crystal wires of cadmium. He found that an oxide coating increased not 
ij only the critical shear stress but also the yield stress of wires that had been 
© considerably extended. The importance of surface effects for glide in metal 
© crystals was also shown by the acceleration of the flow caused by bombardment 
) of the surface with « particles from polonium, which penetrate only a very small 
4 distance into the metal (Andrade 1945). These were the facts known at the 
“| time (1946) when the experiments on surface effects which are here described 
5 were begun,t{ but since then other experiments, cited at the end of this paper, 
© have emphasized the influence of surface conditions. 

The original intention was to find the effect on the mechanical properties 
of a single crystal wire of (a) removing a thin layer of the surface, (b) adding a 
1) thin layer of metal to the surface, both by electrolysis. Experiments on these 
) lines, however, showed that the electrolyte alone, without applied potential 
_ difference, had marked mechanical effects. Accordingly the effect of the contact 
of various electrolytes on the flow of single crystal wires of cadmium was 
investigated: a few experiments have since been carried out with zinc. ‘The 
- general result has been to bring evidence as to the important effect of surface 
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films on the mechanical properties of single crystals, which has also been | 
emphasized by the recent work of Harper and Cottrell (1950). Some electrolytes — 
disrupt a coating of oxide or other surface contamination already existing, and — 
so make the crystal more prone to flow: others lead to the formation of a coating 
on a previously cleaned surface, or thicken the coating on an ordinary surface, — 
in both cases notably raising the critical shear stress. | 


§2. PREPARATION AND TESTING OF SINGLE CRYSTAL WIRES 

The crystals were made from 1mm diameter spectroscopically pure wire, 
drawn by Messrs. Johnson and Matthey from cast rods which had been selected 
as free from cavities. Such cavities, containing air or other gas, remain during 
drawing and lead to uneven surface oxidation in the single crystals. ‘The 
impurity of both cadmium and zinc was less than 0-001% and consisted chiefly 
of lead. A few crystals were made from cadmium of commercial purity; these 
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Fig. 1. Critical shear stress, for spectroscopically pure cadmium single crystal, 
as a function of log rate. 


S(l) =Stress for 1% per hour. 


gave a critical shear stress of 30gwt/mm?, as against 16g wt/mm? for the 

spectroscopically pure metal, both at an initial rate of flow of 0-02% glide per 
minute. ‘This value of the critical shear stress compares well, as determinations 
of critical shear stress go, with the value of 13 g wt/mm? given by Andrade and 
Roscoe (1937) at the same rate of glide. Since critical shear stresses from 
cadmium are quoted at widely different rates of flow, which leads to confusion, 
we give in fig. 1 a curve for spectroscopically pure metal connecting critical ) 
shear stress, expressed in terms of that giving a glide of 1% per hour as unity, | 
with log rate of glide in per cent per hour, calculated from the formula given by | 
Andrade and Roscoe (1937, p. 164), which has been found to give a very | 
satisfactory representation of the experimental facts. At a glide rate of 1% per 
second, for instance, the critical shear stress is 27 g wt/mm? if that for 1% per 
hour is taken as 13 gwt/mm?. The critical stress for the spectroscopically pure 

zinc was found to be 33 g wt/mm? at the rate of glide of 0-02°% per minute. 
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The wires were transformed to single crystals by the travelling furnace 
method of Andrade and Roscoe (1937). The polycrystalline wire was cleaned 
with metal polish and a degreasing agent and then, by means of a special apparatus, 
transferred without admission of air to an evacuated Pyrex tube which had been 
cleaned and baked out under vacuum conditions for half an hour at 400°c. 
In such Pyrex tubes, in which they were held loosely, the wires were sealed off 
at a pressure of 10-*mm Hg. Several of these tubes were contained in a horizontal 
silica tube of about 1 cm internal diameter passing through the furnace, which, 
drawn by a wire controlled by a synchronous motor, travelled along it without 
touching at a rate of about 4cm per hour. The temperature of the furnace was 
maintained constant to within 1:5°c. At either end of the furnace was fixed a 
small compartment, kept cold by running ice water, which ensured that the 
temperature gradient was constant. The whole was carefully shielded from 
draughts. Crystals of any desired orientation were produced by fusing a short 
length of crystal of this orientation to the wire to be transformed. 

Special precautions against strain were taken in the transfer of the very soft 
single crystal from the tube to the extensometer system, which was accomplished 
by means of a specially designed geometric slide. 

The crystal orientation was determined to within 1° by back-reflection x-ray 
photographs and the use of a Greninger chart. In the earlier work the crystals 
were tested at 1 cm intervals along the wire; it was found, however, that crystals 
made in the way described, which gave sharply defined spots in Laue 
photographs, were homogeneous and accordingly in later work only one 
photograph was taken for orientation. 


‘ 


§3. EXPERIMENTS ON PLATING WITH CADMIUM 

The commercial plating solution used* contained cadmium potassium 
cyanide, potassium cyanide, sodium carbonate and caustic soda. In the 
experiments on the effect of plating the surface a control wire was cut from the 
same single crystal wire as the specimen to be plated. ‘This latter was fixed 
vertically in the plating bath on the axis of a cylindrical rotating anode, which 
was provided with vanes for stirring. When the plating process had continued 
for three hours at a current strength that gave a calculated depth of coat of about 
10-5 cm, the wire was removed, washed in distilled water, dried in air, and then 
mounted in the holder with the precautions already described. ‘The crystal 
orientation was measured and the specimen transferred to the extensometer. 

With the control specimen the load was increased by 2 g wt every 5 minutes 
until a slow rate of flow (about 0-01°% to 0-1% per minute) was produced and 
measured. The plated specimen was then tested under the same loading routine 
and the rate of flow compared with that of the plated specimen. In the case of 
eight crystals out of twelve a markedly greater rate for the plated specimen was 
observed, the three highest values of the ratio to that of the control being 
10, 12 and 17. However, the results were irregular and in three cases there 
was actually a decrease of rate, which in the light of subsequent experience 
we attribute to oxidation having taken place in the drying process. Since 
a few control experiments showed that leaving the wires in the electrolyte, 
without passing any current, produced an increase in rate, subsequent 
experiments were carried out under the simpler conditions. 
: * Supplied by Messrs. Canning & Co., Ltd, 
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§4. INFLUENCE OF CONTACT OF ELECTROLYTES 

To avoid subsidiary electrolytic effects the extensometer, shown in fig. 2 
was made of plastic (distrene). The disc AA rested on the edge of a thick glass 
beaker which had been ground smooth, and was held down by weights, not 
shown. To it was attached by three vertical columns, of which one is shown 
at B, the base C, on which was mounted the chuck D, in which the lower end 
of the specimen, about 6cm long, was gripped in two V-shaped cuts by plastic 
screws. The upper end of the specimen was held in a chuck F, attached by two 
threads GG to the pulley wheel P which ran in a ball race, the axis being held 
in the V’s of the support H. The light scale pan S, attached by a single wire 


Fig. 2. Apparatus for straining single crystals in contact with electrolytes. 


to the pulley wheel, counterbalanced the chuck F and carried the weights by 
which the tension was applied. The friction was very small, amounting at most 
to 1% of the total load. ‘The attachment L to the upper chuck carried a fiducial 
mark, allowing the extension of the specimen to be read by a cathetometer 
graduated in mm. 

The method of testing was to increase the load in small steps until a slow rate 
of flow was produced and then, after this rate had been measured for a period of, 
say, 20 minutes, to introduce the electrolyte and measure the subsequent flow. 
The ratio R of the rate after contact of electrolyte to the rate before contact of 
electrolyte was taken as a measure of the effect. ‘This is a very sensitive method, 
as the rate of flow is greatly influenced by small changes in structure or local stress. 

The first experiments were done with the commercial plating solution with 
which the preliminary plating experiments had been carried out. Wires of 
various crystal orientation were used. Immediate increases of rate, with values 
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of R up to 15 were found, but once more the values of R were very variable and 
did not appear to be related to the crystal orientation. It was decided to work 
with simple electrolytes. 

The first electrolyte used was a nearly saturated solution of cadmium 
sulphate. Figure 3 represents a typical experiment, the angle ss, made by the 
crystal axis with the glide plane being 28°.* Thirty-eight experiments were done 
in all, with various values of ys) and with various amounts of pre-glide, from 
0 to 80%. The values of R varied from 1-0 in three cases (which we now attribute 
to the probability that the surfaces were unexpectedly clean) up to 9-6: in 16 cases 
R was between 2 and 3-5. While the effect is marked, and always takes place 
immediately the electrolyte is brought into contact with the wire, it is variable 
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Fig. 3. Flow of cadmium crystal exposed to Fig. 4. Flow of zinc crystal (y= 28°) exposed to zinc 
cadmium sulphate solution, introduced at sulphate solution. The crystal was hardening at 
moment indicated by arrow. Rate of exten- the time, indicated by arrow, when the solution 
sion before introduction of electrolyte, was introduced. 


0-01% per minute: after 0-026 % per minute. 
In this and all related diagrams the arrow 
indicates the moment of introduction of the 
electrolyte. 


in magnitude. However, something can be gathered by taking average values. 
Table 1 shows the average values of R for different ranges of initial angle yo, the 
number of tests from which the average is taken being given. It will be seen that 
there is certainly no significant variation of R with yp, since the slightly higher 
value of R for large values of */, is an average of two tests only. 


Table 1 
ho 0 to 24° 25 to 44° 45 to 64° 65 to 80° 
Average R S30) 2°9 3-0 4-0 
Number of tests 9 20 7 2 


‘The grand average for R is 2-9. 


* It may be remarked that for % greater than 50° local flow may occur, owing to the relation 
between physical hardening and geometrical softening (see Andrade and Roscoe 1937, p. 159). 
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As regards the variation with pre-glide p, table 2 shows average values of R” 
for different ranges of p, from which it seems clear that the effect is smaller for 
small pre-glides, but constant for pre-glides greater than 10%. 


Table 2 4 
Pre-glide (°%) 0 to 9 10 to 24 254049 50 to 80 
Average R 1-4 3338) 3-4 3-4 
Number of tests 8 10 16 4 


It is not surprising that the hardening effect of the surface layer persists after 
glide has taken place, since Roscoe found that it was not very different after 300% 
glide from what it was at very small glide. The significance of the smaller value 
of R for very small pre-glide may be that the layer is more easily attacked by the 
electrolyte if it has already been disturbed by glide. 
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Fig. 5. Flow of cadmium crystal ()=32°) exposed to Fig. 6. Flow of cadmium crystal (%)=48°) exp 
cadmium nitrate solution. Contact of the electro- to zine nitrate solution. ‘There is no in 
lyte produces an initial increase in rate, but after increase in rate: the rate decreases—continuo’ 
five minutes flow has practically ceased. from the moment of contact. 


With zinc crystals a similar increase of rate of flow was produced by the action 
of a concentrated solution of zinc sulphate. Figure 4 is an example, for zinc, 
in which the crystal was hardening at the moment when the solution was 
introduced, so that the value of R is very large. 

A few experiments were carried out with cadmium single crystals and 
cadmium chloride, of the same molecular concentration as the sulphate, as 
electrolyte. ‘They gave similar results. For what it is worth, with results — 
showing a wide spread, the average value of R for the sulphate was 2-9 and for 
the chloride 2:5: the difference is not significant. 

Experiments with cadmium sulphate solution carried out at 50°c showed no 
significant difference from those at 15°c, the average value of R for five trials | 
being 2:1. 

Nitric acid 1 in 500 gave R= 1-6 and nitric acid 1 in 50 gave R=7-0. 

Saturated solutions of cadmium cyanide and iodide had no effect on the flow 
of cadmium crystals, and water, used as a control, likewise had no effect. 
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So far the results recorded with electrolytes have been either negative 
(cyanide and iodide) or have been an immediate increase in rate of flow. With 
cadmium and zinc nitrate solutions, acting on cadmium single crystals, a new 
effect was found. Figure 5 shows a typical result with cadmium nitrate solution. 
There is an immediate increase in rate of flow, as with cadmium sulphate, but 
this is followed within five minutes by an almost complete cessation of flow. 
If the wire is left in the electrolyte for some time and subsequently tested, it is 
found that the load W, needed to produce the rate of flow that was produced 
originally, before the contact of the electrolyte, by a load W,, has considerably 
increased. ‘Thus in one experiment, after 40 minutes immersion, W,/W, was 
1-3, and after 16 hours immersion, in each of two experiments, it was 1-7, which 
is the kind of increase produced in Roscoe’s experiments by an orange coat of 
oxide. In seven experiments on cadmium with cadmium nitrate, 30 minutes 
after the first contact of the electrolyte the rate of flow was reduced to 
one-fiftieth or less, except in one case, where it was one-twentieth. 

The effect of zinc nitrate solution on the flow of cadmium crystals is illustrated 
in fig. 6. The initial increase in rate is absent: instead we have a continuous 
decrease in rate. With this particular wire, after two hours the value of W,/W, 
was 1-4. 

These results, then, are consistent with the view that the Cd ion acts on the 
surface oxide film of an ordinary Cd wire so as to diminish its inhibitory effect 
but that the Zn ion does not do so. The NO, ion builds up a surface layer (oxide 
or hydroxide) which inhibits flow and raises the critical shear stress by some tenths. 


§5. EFFECT OF CLEAN SURFACE 

If, as we had been led to suppose, the effects so far described were attributable 
to surface contamination, it was clearly desirable to carry out some tests with a 
really clean surface. The best way to prepare a clean surface without contact 
with electrolyte is by thermal evaporation and accordingly cadmium single 
crystal wires were cleaned by heating zm vacuo (pressure <10-*mm Hg) in the 
apparatus described by Andrade and Randall (1950). A layer about 4x 10-?cm 
was removed, leaving the surface covered with pits, as described in the paper just 
quoted. 

In a typical experiment a length of 6-5 cm was removed from a single crystal 
some 30cm long and the remainder was surface-cleaned by evaporation. With 
the length removed a normal extension experiment was carried out with the 
apparatus of fig. 2, cadmium sulphate being used as the electrolyte. A section 
of the cleaned crystal was then placed in the extension apparatus and the load 
adjusted to give the same initial rate as with the uncleaned specimen, the interval 
between removal from the vacuum and the test being kept as small as possible. 
After 20 minutes flow at the chosen small rate the liquid was introduced. 

Figure 7 shows a typical result. The rate ratio R is 1-5 for the normal surface: 
with the cleaned surface the introduction of the electrolyte does not affect the 
rate. In the following table are summarized the results of the eight experiments 
carried out. It will be seen that in no case was there any change of rate with the 
clean surfaces, while with the uncleaned crystals there was in all cases an increase 
in rate. That the rate ratio never exceeded 1-5 may be taken as an indication 
that the surface contamination was never bad. 
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In the case of the first two specimens there was a pre-glide of 30%, with the 
other specimens no pre-glide. | 

io, Xp are the angles initially made by the wire axis with glide plane and glide 
direction; R is the ratio defined on p. 448 for the ordinary crystals and R, for the 
crystals surface-cleaned in vacuo; o and o, are the respective critical shear stresses. 
pertaining to the same initial rate of glide. 
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Fig. 7. Behaviour of surface-clean and normal cadmium crystal (6,= 22°) contrasted. 


Table 3 : 
to Xo R Re 0,/o 
20 20 1-4 1-0 1-8 
26 26 1-1 1-0 sh? 
10 10 152 1-0 15 
10 12 ila 1-0 isp2 
ay B35) 1°5 1:0 i133) 
33 33 1-4 1-0 1:3 
40 46 1-4 1-0 Lee 
42 42 1-4 1-0 1:2 


It will be observed that the critical shear stress is always higher for the 
clean-surface wire than for the normal wire. Andrade and Randall (1950) have 
attributed the pits formed by thermal etching to local spots of high free energy, 
where the vapour pressure is abnormally high. An alternative statement would 
be that thermal etching results in the removal of atoms which, owing to local 
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irregularity, are comparatively loosely bound. Since such local irregularity is 
likely to act as a source of dislocations, the rise of critical shear stress is 
understandable: local sources of weakness are boiled away. 


§6. POLYCRYSTALLINE WIRES 

A few experiments were made to see if there was any effect with polycrystalline 
wires. Tests on commercial cadmium wire of 1mm diameter, with about 
40 crystal grains to the diameter, failed to reveal any influence of the plating 
solution which had so marked an effect with the single crystal wire. The flow 
of polycrystalline zinc wire was likewise uninfluenced by zinc sulphate. 
Spectroscopically pure cadmium wire was annealed at 300°c until the grains had 
so grown that there were only one or two to the diameter. With this wire, owing 
to the large grain size, the extension was somewhat irregular; however, neither 
plating solution nor cadmium sulphate solution produced any observable 
increase in rate. 

§7. DISCUSSION 

The general result of these experiments, of which a preliminary note was 
published some time ago (1948), is to emphasize the great influence of surface 
layers of oxide, and/or possibly other salts, on the mechanical properties of 
single crystals, an effect which is particularly striking with silver crystals, as 
reported by Andrade and Henderson (1951). The absence of effect with a 
really clean surface, which has been confirmed by Menter and Hall (1950), seems 
to render it highly probable that the increased flow produced by suitable 
electrolytes is due to disruption of a film which tends to inhibit the propagation 
of sensitive dislocations from the surface, while the hardening action of both 
zinc nitrate and cadmium nitrate, which takes some time to establish itself in 
both cases, is due to the formation of a thin surface layer, which, according to 
Menter and Hall (1950), consists of cadmium hydroxide. 

The immediate increase of flow caused by a solution of cadmium nitrate, 
which is not followed until later by the hardening, points to the cadmium ion 
as the effective agent in disrupting the surface film to which resistance is due. 
The contrast with the action of zinc nitrate, which produces from the start a 
progressively increasing hardening, points to the NO, ion as the effective agent 
in producing the hydroxide film. 

It appears, then, that the action of solutions of cadmium sulphate and chloride 
is due to the cadmium ion. We are not fully clear as to the mechanism of the 
action, but probably the cadmium ion combines with the oxide or hydroxide 
component of the surface layer and so removes it into solution. 

The lack of action of concentrated solution of cadmium iodide and cyanide 
can be explained, on the basis of views here put forward by the well-known fact 
that concentrated solutions of these salts give rise to very few Cd?* ions. As 
shown by Riley and Gallafent (1932) the bulk of the metal in such concentrated 
solutions is present either as the undissociated molecule or as CdI* ions. 
A further point is that the greater ionization of the sulphate and chloride leads 


to greater hydrolysis 
CdSO,+2H,O = Cd(OH), + H,SO, 


so that they would have a greater solvent action on a film of cadmium oxide than 
the iodide and cyanide. 
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Work by Raether (1950) on zinc crystals with freshly cleaved surfaces shows _ 
that any such surface exposed to air is covered by an oxide layer some 20-404 — 
thick and indicates that electrolytically polished surfaces also become 
contaminated. Roscoe’s work shows that a layer of 20-404 thickness is sufficient 
to raise appreciably the critical shear stress of cadmium single crystals. The 
strength of a crystal may be permanently lowered by electrolytic action, as shown 
by our preliminary results on cadmium and by the work of Harper and Cottrell 
(1950) on electropolished crystals of zinc. This reduction in mechanical resistance 
may, perhaps, be attributed to a greater ease with which a travelling dislocation 
can pass out of a crystal at a perfectly clean surface, in which case it 1s 
understandable that further immersion in electrolyte would produce no effect, 
despite the presence of a film some 20-404 thick, in accordance with the results 
of Harper and Cottrell. 

The lack of effect with polycrystalline metal is easily understood on the 
views that we have put forward, since only the surface grains would be affected, 
but not on the views of Rehbinder and Wenstrém (1937), who consider that a 
deep penetration of the active agent takes place. Grain boundaries would seem 
to offer great facilities for such action. 
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Ferromagnetic Resonance in Manganese Antimonide 


Ferromagnetic resonance absorption has been studied in powdered specimens of MnSb 
at a wavelength of 1:26 cm and in the temperature range from 20°c to the Curie point | 
(315°c). ‘The absorption measurements were performed as previously indicated (Adam | 
and Standley 1950), the powdered specimen under examination being formed into a thin 
plate, using a suitable low-loss adhesive. ‘The applied steady magnetic field, H,™®*, for 
maximum power absorption was determined at each temperature. 

Kittel (1948) has shown that the resonance condition is given by 
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Here v is the frequency of the microwave field H,,, g is the spectroscopic splitting factor, 
N,, Ny, N; are the demagnetizing factors of the specimen in the three coordinate directions, 
and the other symbols have their usual significance. (N,,—N,)I was found to be negligible in 
these experiments, and an oscillation magnetometer of the type described by Griffiths and 
MacDonald (1951) permitted the measurement of the factor (N,—N,)I throughout the 
required temperature range on the actual specimens used in the absorption measurements. 

Typical results are shown in the table, and it will be noted that, above about 100°c, as 
the temperature rises both the g-value determined from the above equation and the half- 
width of the absorption curve decrease. Five different samples all showed the same behaviour, 
with only small deviations from the mean values given in the table. 


Temperature (°c) 20 65 Oh; 140" 1 7200 200 = 2205 250 00 
Half-width (Oe) 5000 5500 5700 5600 5000 4000 3300 2500 1100 
g AN Meal oe Hi ead 1 aur Feeney elas let RN aie oF 


Griffiths (1951) and MacDonald (1951) have shown that anisotropic strains in the 
specimen modify the above resonance formula, but that the form of the equation remains 
unaltered if the terms which account for shape anisotropy are written CN Nai) f anid 
(N,, —N,/)J and are taken now to include also strain anisotropy. The oscillation magneto- 
meter yields the modified term (N,,’— N,’)J, and there is therefore no obvious experimental 
condition which would account for the observed change in g-value with temperature. 

The above results are at variance with those found by Bloembergen (1950) for nickel and 
supermalloy. In these materials the g-value is unaffected, within experimental error, by 
temperature changes and the absorption curves broaden appreciably as the temperature is 
raised. A recently published note by Okamura, Torizuka and Kojima (1951) indicates that 
results, somewhat similar to those reported here, have been found in cobalt ferrite. 

Kittel (1949) and Van Vleck (1950) have shown that spin-orbit interactions can lead to a 
spectroscopic splitting factor greater than the spin-only value of 2:00. However, in the 
temperature range from 20°c to 300°c the change in lattice configuration should be small, 
so that the contribution of the orbital motion to the g-value should not appreciably alter, and 
hence a constant g-value would be expected. The change in line width with temperature in 
MnSb and in cobalt ferrite could be explained in terms of a decrease in the anisotropy 
constant K with increasing temperature. K is not known for MnSb, but if it is high at room 
temperature (as in MnBi and cobalt ferrite) then the observed decrease in g-value follows the 
probable decrease of the anisotropy forces. Although the experimental accuracy diminishes 
near the Curie point, due to reduced absorption, g appears to approach the value 2-1 rather 
than 2:0. Here the anisotropy forces are expected to be small, and the departure from the 
-spin-only value is similar to that found in such materials as supermalloy and nickel. 


Department of Physics, G. D. Apam. 
University of Nottingham. K. J. STANDLEY. 
23rd March 1952. 
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Low Remanence and the Temperature Variation of 
Permeability of Silicon Iron Alloys 


According to simple domain theory the remanent intensity of magnetization IR of a 
polycrystalline ferromagnetic should have a value equal to 0:5 Ig when the material 
possesses uniaxial anisotropy, and 0-83 Ig or 0°87 Ig according as the direction of easy 
magnetization is the [100] or the [111] direction respectively. In most materials these 
predictions are fairly well supported by experimental observations, but there are some 
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which possess much lower remanence than that given above. Bozorth (1947) has suggested 
three possible mechanisms whereby these low values may be brought about, and he 
considers that the most likely explanation is that in magnetically soft materials the domain 
wall energy is sufficiently small to allow the substance to demagnetize itself by the formation 
of new domain boundaries as soon as the applied field is removed, thus giving rise to quite 
large regions inside the material where the flux is closed to a considerable extent. The 
greater the extent of this flux closure, the smaller the remanence becomes. 
The expression quoted by Bozorth for the wall energy per cm? is 


3 K\12 (ro\1!2 
W=2 (5 Ao+ 5) (=) 


where 2 is the saturation magnetostriction, o the stress acting on the material, K the 
magnetocrystalline anisotropy coefficient, k Boltzmann’s constant, @ the Curie temperature 
in degrees Kelvin and a the interatomic distance. From this equation it is evident that 
W decreases with increasing temperature. Now it is possible that flux closure may occur 
to a lesser extent even in the presence of a small applied field, and if this mechanism occurs 
at any stage in the course of magnetization, then at a given field strength the intensity of 
magnetization will vary with temperature in a manner different from that to be expected 
when flux closure effects are absent. In particular, if the process occurs to any great extent, 
the permeability may decrease with temperature in contrast to the behaviour usually 
observed. 

It is to be expected that the materials which would exhibit this effect in the greatest 
degree would be those for which the saturation intensity, and hence the self-demagnetizing — 
energy, is large, and for which the crystal anisotropy is small. These properties are 
possessed by the iron-rich Si—Fe alloys. 

Measurements made in this laboratory have shown that for 4°% Si-Fe the remanence 
is about 400 gauss, so that Ip/Ig=0-25. According to Bozorth’s figures there are 
indications that the remanence decreases with increasing silicon content, the ratio Jp/Ig 
for a 65% Si-Fe alloy being only 0:11. Now it is well known that these alloys exhibit an 
anomalous temperature variation of permeability, the permeability decreasing with 
temperature at low field strengths in contrast to the behaviour of pure iron in this respect. 
Spooner (1926) has shown that the addition of silicon to iron causes du/dT to decrease 
and to acquire negative values at certain inductions, du/dT becoming more negative 
the greater the silicon content. For example, at room temperature and for B=10000 gauss 
_du/dT in gauss per oersted per °c is —0-00012, —0-0008, —0-0018 for silicon—iron 
containing 0:9% Si, 2-294 Si, and 4% Si respectively. There appear to be no similar 
measurements for 6:5°%, Si-Fe, but observations by Zaimovsky (1941) indicate that for 
12-14% Si-Fe du/dT is also negative in low fields, and the same worker has found similar 
behaviour in iron-rich Fe—Si—Al alloys. 

These observations suggest that low remanence and negative temperature coefficients 
of permeability are closely connected and that both can be explained in terms of internal 
flux closure. The bearing of such a magnetic process on the low field magnetocaloric 
effect, which would give rise to anomalous temperature changes in low fields of the same 
sign as the high field effect, has been discussed by Stoner and Rhodes (1949) and has 
recently been confirmed by Bates and Marshall in unpublished experiments on Si-Fe 
alloys to be described shortly. 

I am indebted to Professor L. F. Bates for his interest and advice and for allowing me 
to quote the above results prior to publication. 


Department of Physics, E. W. Lee. 
University of Nottingham. 
19th March 1952. 
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Strain Patterns in Toughened Glass 


When plane polarized light is transmitted by a sheet of toughened glass at large angles of 
incidence, the strain pattern may be clearly seen without the usual analysing screen or prism. 
The effect may be observed simply by looking through an inclined sheet of strained glass 
towards clear sky, and is sometimes noticeable in motor car windscreens. 

The Fresnel formulae show that approximately 25°/, of the light transmitted by a single 
glass surface is plane polarized for an angle of incidence of 80°. A plate of toughened glass 
may therefore be considered to contain its own analyser, the glass—air surface, if it is 
sufficiently inclined to the incident polarized light. An auxiliary, strain-free plate, held 
parallel to the toughened panel, should augment the analysing effect and increase the 
contrast under which the pattern may be seen. ‘The two added surfaces would then 
contribute a degree of polarization intermediate between that found from the Fresnel 
formulae and the corrected value allowing for internal reflection (Gaviola and Pringsheim 
1924, Tuckerman 1947). The corresponding increase in contrast should be roughly twice. 
Polarization at the front surface of the toughened plate is insufficient, by itself, to lead to the 
formation of a strain pattern. 

Experiments have shown that sufficient contrast is obtained with angles of incidence 
greater than about 70° for the pattern of a strained glass specimen to be seen in plane polarized 
light, without the usual analyser. The addition of an auxiliary parallel glass plate has led 
to a marked increase of the pattern contrast. By rotating the plane of polarization relative 
to the inclined panel the intensity of illumination has been found to show two maxima and 
two minima per revolution, at each point of the pattern. This is in accordance with the 
formula, given by Schuster, cos? (8—a)—sin 2a sin 28 sin? 6/2, where a, 6 are the angles 
between the optic axis and the directions in the polarizer and analyser respectively and 6 
is the phase difference. The formula is here modified for the case B=constant, and 
for points of equal phase difference between components the intensity varies as 
J+ cos 2a+/ sin 2x, where j=4 and k and / are constants which depend on the values 
chosen for 8 and 6. With incident white light the pattern is coloured in the usual manner. 

Field tests have been made with the light reflected from large surfaces such as roads or 
runways, with skylight at altitudes exceeding 20 000 ft. and with light reflected from layers 
of light cloud or haze below an aircraft. All tests have shown the importance of the 
strain pattern effect in oblique windscreens incorporating toughened panels. Very pro- 
nounced patterns were observed against layers of thin haze and this effect was still very well 
marked during twilight. Since analysis occurs at the glass—air surface the effect can only be 
obviated by lessening the strain pattern in the toughened glass. 

Thanks are due to the Chief Scientist, Ministry of Supply, for permission to publish this 
note. ‘ 


Royal Aircraft Establishment, J. M. Natsu. 
South Farnborough, Hants. E. R. WEBB. 
3rd April 1952. 
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On the Sampling of Water Droplets in Natural Clouds 
and in Radiation Fogs 


In their recent review of the microphysics of clouds, Mason and Ludlam (1951) 
include a section on the determination of cloud droplet size which calls for certain 
comments. It may be that some facts escaped the notice of the authors or that they were 
omitted as being unimportant, but their statement that “ a glass slide coated with vaseline 
or magnesium oxide may be exposed to the cloud elements ” suggests an over-simplification 
of the problem of size determination. 
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The mixture in which water particles are trapped and entirely confined must fulfil 
certain requirements, for example its density should be nearly the same as that of water, 
the matrix must be smooth and highly viscous (Benedicks and Sederholm 1944) but of 
refractive index different from that of water to ensure good contrast. Pure vaseline, as . 
proposed by the authors, is definitely not suitable as a matrix for slide-surfacing for the 
following reasons, clearly stated by Houghton and Radford (1938) : the surface of the — 
vaseline layer on the glass is not smooth, and when viewed under dark ground illumination — 
it is streaky and irregular, so that some of the captured drops are always distorted and the 
contrast is very poor. Further, the density of vaseline is 0-873 g/cm® at 23°c. 

The problem of proper surfacing was solved quite satisfactorily by Fuchs and 
Petrjanoff (1937) who found that good results can be obtained by collecting the droplets 
on glass slides coated with a freshly melted mixture of light mineral oil and vaseline. The 
proportion of the ingredients has to be varied according to air temperature and is 
3:1 at 20°c. — It is quite clear that the light mineral oil is the more important ingredient— 
its refraction is higher than that of the water. 

This method was first used in this country by the present writer during his war-time 
investigations on natural clouds and radiation fogs (Mazur 1943; see also Kraus and 
Smith 1949, Perrie 1950 and Best 1951). The light mineral oil was previously saturated 
with distilled water and every effort made to obtain a bubble-free film mixture. In order 
to avoid a piling-up of droplets into layers when being captured on the plate, the duration — 
of sampling has to be very short. The instrument used was known as the aircraft impactor 
and consisted essentially of a high-speed shutter with the glass plate covered either with 
a uniform film of the Fuchs—Petrjanoff mixture or, alternatively, with a magnesium oxide 
layer on which the impinging droplets made craters. / 

The water droplets collected within the oil—-vaseline mixture can be preserved for 
many hours without showing any perceptible diminution in size or change of shape. 

The method adopted for the collecting of drops in natural clouds was to coat the cover 
glasses with the fresh mixture just before the aircraft take-off and then to warm them. 
before use in flight by putting them on to the flat surface of a lagged vessel filled with hot 
water at about 45—50°c. 

Some authors have obtained satisfactory results using other mixtures: (1) Hauser and 
Strobl (1924) used neat’s-foot oil and glycerine; (11) Diem (1942) used two layers of different 
oils of different densities, collecting the drops on the interface between the layers; 
(iii) Bricard (1939) used mineral oil to which organic acids had been added in order to 
reduce the surface tension to render easier penetration. 

The magnesium oxide method is only a modification of the method developed by 
Strazhevsky (1937), who used a uniform layer of soot smoked on to a slide during his 
investigation on the atomization of liquid fuel. ‘The magnesium oxide layer is very fragile 
and can be damaged by high-velocity air jets during exposure to the droplets. ‘The author 
found that the oxide layer is blown off the plate at a speed of about 200 m.p.h. The 
method is of little value for droplets smaller than 10 wu. 

Incidentally Pauthenier and Brun (1940), mentioned by the authors, have also used 
the well-known ‘stain’ or ‘absorbent paper’ method. 


British Rayon Research Association, J. Mazur. 
Manchester. 
9th November 1951. 
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The Dispersion of Electron Beams in Gases 


When a beam of electrons passes through a gas it is scattered mainly elastically, because 
self-repulsion (Watson 1927) is negligible. The size of a spot produced by the diffuse beam 
on a fluorescent screen can be measured by a radius 7; (in cm) which corresponds to a 
point on the screen having 1/e of the intensity at its centre (Lenard 1927); theoretically 

ny Sel hd BAI e ey a Se eee (1) 
Here ¢ (of order 10° for N,) depends on the gas, p is the pressure in mm Hg, V,, the 
accelerating potential in volts and d the distance travelled in cm. For p>10-? and 
V,=10* ...10* a spot should not form. 

It was found very early that ionization by an electron beam in a rarefied gas (e.g. in 
He or A) concentrates the beam (Arnot 1930, 1931, Cosslett 1950). A treatment of this 
* gas-focusing ’ was not very successful and the general results could not be compared 
with experiment (Scherzer 1933). In all these cases (Briiche 1932) the gas pressure used 
never exceeded 10-*mm Hg and this has been accepted as a practical upper limit. 
Fortunately this conclusion was unjustified. 

We have discovered that an electron beam accelerated to 3-10 kv in nitrogen remains 
well concentrated even at pressures of about 1mm Hg. We believe that 1 mm Hg does not 
represent an upper limit and that the nature of the gas is not critical. 

The ultimate cause of the concentration of the beam is the presence of positive ions 
in and around the beam. The fast electrons produce at this pressure a small number of 
ion pairs along the path of the beam.’ The slow secondary electrons and positive ions so 
formed diffuse radially outwards. ‘The electrons, being originally faster than the ions, 
set up a distribution of concentration such that at the axis of the beam the concentration 
of ions slightly exceeds that of the electrons; as a result a positive net space charge develops. 
Since secondary electrons and ions are produced in equal numbers, a negative net charge 
develops at larger radii. At still greater radii the relative velocities of electrons and ions 
become so small that the particles recombine readily in the gas and their concentration 
decreases more rapidly. As the space-charge distribution develops the radial velocities 
of electrons become smaller and those of the ions larger. In equilibrium their speeds at a 
given radius are equal and so are the numbers crossing the corresponding cylindrical 
surface. In the following arguments the effect of the negative space charge due to the 
fast electrons is neglected. 

We have made some rough calculations in order to test whether this picture gives 
consistent results. "The condition to be satisfied is that the net positive charge in the beam 
is large enough to restrain the fast electrons from being scattered out of the beam. In 
addition we have to assume that the average random velocity of the secondary electrons 
is larger than the corresponding velocity of the fast electrons. 

An electron of energy 5 kev produces about 1 ion pair per cm at 1 mm Hg (von Engel 
and Steenbeck 1932). Thus with a beam current of 1 microampere, 6 x 10!” ion pairs are 
produced per sec per cm length of the beam. If the loss by ambipolar diffusion is balanced 
against this rate of production the equilibrium concentration | found is of order 
1014 ion pairs/em*. A similar value is obtained if the concentration is determined by 
balancing recombination against diffusion, using a recently measured value for the 
coefficient of recombination (Biondi and Brown 1949). The above value of the 
concentration is large enough to justify the assumption of ambipolar diffusion (Loeb 1939). 
At the same time it is not excessive and thus the concentrations of positive ions can be 
slightly different from that of the electrons. This is particularly likely to occur at the 
core of the beam. 

If ambipolar diffusion operates, the net positive charge should be only a fraction of the 
total charge present. We have calculated therefore the space charge needed by assigning 
to the primary beam an average energy of random motion of order 10 ev, a value which 
can be obtained from the classical treatment of elastic scattering. Since the beam has a 
radius of less than 1 mm the potential at this radius must be about 10 v negative with 
respect to the axis if the fast electrons are to be retained within a beam of ee radius. 
The space charge density which gives this potential is of order 10''ions/cem*. These 
results appear to be consistent and reasonable. 
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The axial motion charges depend on the potential of the isolated screen with respect 
to the anode. In high vacuum this potential may be large (McKeag and Welch 1951, 
McKay 1948) but in a gas the screen potential must be so small that the radial charge 
distribution is not disturbed. 

Gas focusing at low pressure is a fundamentally different phenomenon. The radial 
motion of the secondaries occurs essentially in a vacuum and,so the walls become important 
(Arnot 1930). 

The permission is acknowledged of the Director of the Electrical Research Association 
to publish this note. Professor Lord Cherwell has kindly placed the facilities of this 
laboratory at our disposal. 


Clarendon Laboratory, PR. E. orrerpes 
Oxford University. A. VON ENGEL. 
28th March 1952. 
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Microwave Faraday Rotation in Liquid Oxygen 


Observations of the microwave Faraday effect have been reported by various workers 
(see references) on paramagnetic and ferromagnetic materials, and theoretical discussions of 
the ferromagnetic case have been given by the present writer (1951) and by Hogan (1952). 
Oxygen gas and liquid are well known to be paramagnetic, and it has been of interest in 
connection with rotation measurements on other substances at liquid oxygen temperature 
to know what contribution might be made to the observations by atmospheric oxygen in the 
waveguide. ‘The expression given previously (Roberts 1951) for the rotation in ferromag- 
netics may be applied to paramagnetics within the same approximations provided that we 
use the appropriate value for the intensity of magnetization (in the applied field). Thus ina 
field of 1000 oersteds we estimate for liquid oxygen just below its boiling point at atmospheric 
pressure (assuming a dielectric permittivity of about 1-5) a rotation of about 0-6 deg/cm. 
Measurements in these Laboratories with sample lengths of about 15 cm at a frequency of 
10 kMc/s have failed to show any rotation as great as 0:2°% of this estimate, either for the 
liquid or for gaseous O, just above its boiling point. 

The almost complete suppression of the Faraday effect in this case may be attributed to 
the very strong magnetic damping of the spin precession caused by the molecular collisions. 
Such collisions will be exceptionally effective in the case of the oxygen molecule in the micro- 
wave frequency range, on account of the inherent (magnetic) molecular resonance which 
occurs at about 60 kMc/s. 


Post Cffice Research Station, F. F. RoBErtTs. 
Dollis Hill, London. 
24th April 1952. 
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The Influence of Heat Treatment on Magnetic Viscosity 
in Permanent Magnet Alloys 


In the development of the magnetic properties of dispersion hardened permanent magnet 
alloys, the usual form of heat treatment consists of quenching the material from about 
1 300° c in air or water and then maintaining it at a lower annealing temperature of about 
650° c for a period of 2 to 4 hours (Bates 1951). The alloys are stable in a single phase 
_ System at 1300° c, but have a transition temperature, which lies between 700° and 900° c, 
below which a two-phase system is the most stable form. Thus in the case of Fe,NiAl, 
which may be considered as the prototype of many permanent magnet alloys, Bradley and 
Taylor (1938) have shown that a single body centred phase of Fe,NiAl is stable above 900° c. 
Below 800° c the alloy tends to form a two-phase system, the component phases both being 
body centred and having the approximate compositions FeNiAl and Fe. In the normal 
heat treatment of these alloys, the two phases have insufficient time to separate completely 
but occur in the material as islands or as small inclusions. It is the presence of the inclusions 
which most probably accounts for the high coercivities of the alloys (Néel 1946). 

The metallurgical processes which occur during this type of heat treatment can be 
envisaged in terms of the two mechanisms reviewed by Fisher, Hollomon and Turnbull 
(1948), viz. (a) nucleation, that is, the formation of nuclei in the parent matrix from which the 
inclusions may eventually develop and (6) the growth of nuclei into islands and inclusions. 
The considerations of Fisher et al. show that to a first approximation, during a sufficiently 
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rapid quench from 1300° c, only process (a) occurs to any significant extent. During the 
subsequent anneal the conditions are such as to inhibit the formation of new nuclei but the 
growth of existing nuclei, process (5), is favoured. ‘To a good approximation, therefore, the 
number of inclusions and their size can be independently controlled by the type of quench 
and time of anneal respectively. 

The coercive force of a specimen is a function of the product of the size and number of 
the inclusions (Néel 1946) and consequently, it is to be expected that widely differing types 
of heat treatment can lead to the production of specimens of the same coercive force. This 
behaviour is observed with alnico, e.g. the same coercive force can be achieved either by a 
relatively slow quench from 1250° c, or by a rapid quench in cold water followed by an 
annealing period at 700° c. On the other hand, the phenomenon of magnetic viscosity in an 
alnico specimen is a unique function of the previous heat treatment. ‘The magnetic viscosity 
may be described in terms of the parameter S' defined in a previous paper (Street and Woolley 
1949) in which it is shown that the time variation of intensity of magnetization AJ of alnico 
specimens is of the form AJ—S In t+ constant. In this expression, ¢ is the time which has 
elapsed after a sudden change in the magnetic field applied to the specimen has been made. 

The figure shows the results of measurements of S and H, on alnico specimens which were 
subjected to two different types of heat treatment: (a) cooled from 1250° c at different 
rates, varying from water quenching to cooling the specimen at the natural rate of cooling of 
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the furnace with the current switched off; (5) cold water quench followed by successive 
anneal at 700° c for various total times varying from 5 minutes to 140 hours. H 

Work on the lines indicated above is proceeding with a view to explaining the observed 
phenomena in terms of the nucleation theory of precipitation from solid solution. It is 
hoped to publish an account of this work at a later date. 

A further point arises from the results given in the figure, which extend the rather 
similar measurements which have been briefly reported by Bulgakov and Kondorskii (1949). 
In the latter work, it is concluded that the smaller the coercive force of alnico the greater is 
the magnetic viscosity. The present results show that this statement is only true for values — 
of coercive force exceeding about 150 oersteds. From the figure it can be seen that speci- 
mens of coercive force smaller than 150 oersteds have, in fact, smaller magnetic viscosity the 
smaller the coercive force. 


Department of Physics, R. STREET. 
University of Nottingham. J. C. Woo.LLey. 
15th April 1952. P. B. SmiTH. 
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ABSTRACTS FOR SECTION A 


The Thermal Conductivity of Metals at Low Temperatures: I—The Elements 
of Groups 1, 2 and 3, by K. MENDELSSOHN and H. M. RosENBERG. 


ABSTRACT. 'The thermal conductivity of Cu, Ag, Au, Mg, Cd, Al, In and one poly- 
crystalline and two single crystal specimens of Zn have been measured in the range 2—40° k 
Each specimen shows a maximum in the conductivity in this range. The values of the 
coefficients « and f in the equation for the thermal resistance R=a7T?+{/T have been 
calculated for each specimen. 


The Thermal Conductivity of Metals at Low Temperatures: II—The Transition 
Elements, by K. MENDELSSOHN and H. M. RosENBERG. 


ABSTRACT. ‘The thermal conductivity in the range 2-40°K has been measured for 'Ti, 
Mn, Fe, Ni, Zr, Cb, Mo, Rh, Pd, Ta, W, Ir, Pt, U, and for a Pb single crystal. As a rule, 
the metals of group 8 ae a much higher conductivity than have the metals of the earlier 
groups. The values of « and f in the equation for the thermal resistance R=aT?-+- f/T have 
been calculated for each specimen. The experiments on Cb indicate that the normal state 
was not reached in a field of 3300 gauss at temperatures below 5°K. 


On the Limiting Vibrational Frequencies of an Ionic Lattice, by M. BLACKMAN. 


ABSTRACT. The vibrational frequencies of an orthorhombic ionic lattice, of the 
thallium fluoride type, have been worked out in detail for wavelengths long compared with 
the crystal spacing though short in relation to the crystal size. It is shown that these 
frequencies are not unique but depend in general on the direction of the wave in the crystal— 
in contrast to what happens in cubic crystals. The effect of this behaviour on the vibra- 
tional spectrum is discussed briefly. 


Free Energy of the Double Layers of Two Parallel Plates in a 1-2 Electrolyte, 
by S. Levine and A. SupDDaBY. 


ABSTRACT. By making use of the Poisson—Boltzmann equation an expression is obtained 
for the free energy of the electric double layers of two parallel plates immersed in a large 
volume of an electrolyte of the 1-2 valency type. ‘Two cases arise depending on whether 
or not the divalent ion and the surface charge have the same sign. ‘The formula for the 
free energy is expressed in terms of elliptic integrals. An approximate form for the free 
energy, which is more convenient for numerical computations, is also obtained. This is 
particularly suitable at large surface potentials and is derived for a symmetrical electrolyte as 
well as for the asymmetrical ones. Numerical tables of the free energy as a function of 
plate separation are given and some applications are briefly indicated. 


The Near Ultra-Violet Band System (D-X) of Silicon Monosulphide, by 
A. Lacerovist, G. NILHEDEN and R. F. Barrow. 


ABSTRACT. The rotational analysis of the near ultra-violet band-system, p!IJ—x!X*, has 
been extended. ‘Twelve bands have been analysed, namely, the 7,1, 6,0, 6,1, 5,0, 
5,1, 4,0, 4,1, 3,0, 2,0, 1,0, 1,1 and 0,0 bands. The constants derived for p1#II are 
B, =0-2664,—0-0021, (v’ +4); Dy’ =[2:91—0-03(0’+ $)] x 10~, 
G, =513-1,(v’ +4)—2°95(v' +4)"; ve=35026°8, cm™}; re =2:058, A. 
A number of perturbations have been found in the upper state. ‘They appear to arise from 
interactions with two electronic states, a °X” state and a singlet state, Il or 1A. It has 
-not been possible to determine the absolute numbering of the vibrational levels of the 
perturbing states, but in all other respects fairly complete information has been obtained 
about the 3=- state. ‘The perturbations are similar to those already known in the spectrum 
of CO, and the analogy between the electronic states of the molecules is thus 


strengthened. 
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Dissociation of Molecular Hydrogen Ions in the Cyclotron, by K. E. A. Errat. 


ABSTRACT. The dissociation of H,* ions by collision with gas molecules greatly reduces 
the current of these ions that can be obtained from a cyclotron. The current has been 
measured at two distances from the ion source corresponding to energies of 9 and 18 Mev 
and with various pressures of air and argon in the tank; the loss of ions in expanding from 
the smaller to the larger radius is somewhat less than is predicted by Salpeter’s calculations 
of the dissociation cross section. Assuming this cross section to be inversely proportional 
to the kinetic energy of the ion, the product cH is found experimentally to be 
(1:9+ 0:4) x 10-16 cm? Mev per ‘atom’ of air and (7:0+1-4) x10-' cm? Mev per atom 
of argon. 


The Increase of Ionization with Momentum for Energetic Cosmic-Ray Particles, by 
J. Becker, P. CHanson, E. NaGerotte, P. TREILLE, B. T. PrRIcE and ~ 
P. ROTHWELL. 

ABSTRACT. An experiment is described in which the momentum and specific ionization 

of cosmic-ray mesons were simultaneously measured, the momentum by means of a cloud 

chamber in a magnetic field, and the specific ionization by rectangular proportional counters. 

The difficulties inherent in the use of proportional counters for such measurements are 


discussed. 
The experimental results are in satisfactory agreement with theoretical predictions of an 


increase of specific ionization for very energetic “-mesons. 
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